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FOREWORD 
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ABSTRACT 


Five  gas  turbine  engine  mainshaft  roller  bearing  configurations  were 
investigated  for  capability  of  sustained  performance  at  DN  values 
(Bore  in  mm  x  Speed  in  RPM)  f  *om  2  million  to  3.5  million  and  normal  operat¬ 
ing  temperatures  to  600° F.  A  unique  Dual  Diameter  Roller  was  selected 
for  the  final  analysis  of  stress  and  lubrication  parameters,  design  and 
fabrication.  A  140  mm  Bore  Dual  Diameter  Roller  Bearing  operated  success¬ 
fully  for  30  min.  continuous  operation  at  25,000  RPM  (3.5  million  Dtl) 
with  stabilized  outer  race  temperatures  above  525°F.  Lubrication  was 
with  Polyphenyl  Ether  5P4E  in  an  air  environment. 
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SECTION  I 


INTRODUCTION 

The  objective  of  this  research  program  was  to  develop  a  long  life, 
high  speed  roller  bearing  design  for  future  high  performance  turbine 
engines.  The  bearing  must  be  capable  of  sustained  performance  at  DN 
values  from  2.0  x  10°  to  above  3.5  x  10  and  normal  operating  tempera¬ 
tures  to  60C  F  while  maintaining  the  rotating  shafts  within  the  neces¬ 
sary  tolerance  for  optimum  compressor'  and  turbine  performance. 

The  trends  in  turbine  engine  technology  are  toward  increased  rota¬ 
tional  speed,  larger  bearing  diameters  and  increased  bearing  temperature. 
It  is  anticipated  that  bearing  DN  requirements  (bearing  bore  mm  x  rpm) 
will  grow  from  the  current  2  x  10^  or  less  to  greater  than  3.5  x  10^  in 
the  1980's.  Efforts  to  increase  bearing  DN  capability  are  currently 
aimed  at  ball  bearings.  Approaches  being  investigated  for  ball  bearings 
include  hollow  balls  and  compliant  races.  The  roller  bearing  has  re¬ 
ceived  very  little  attention.  Efforts  which  have  focused  on  roller  beat¬ 
ings  indicate  that  variations  in  roller  bearing  geometry  such  as  the 
dual  diameter  roller  and  hollow  rollers,  show  promise  for  significant 
increases  in  bearing  performance.  Improved  ball  and  roller  bearing 
capability  is  required  for  future  systems. 

Known  current  roller  bearing  types  were  investigated  in  order  to 
select  a  preliminary  configuration  for  further  development.  The  selec¬ 
tion  was  based  on  the  potential  of  each  type  to  meet  the  capabilities 
shown  in  Table  I.  A  bearing  bore  of  140  mm  was  selected  for  this 
preliminary  configuration  study. 

The  results  of  the  preliminary  configuration  study  resulted  in  the 
selection  of  the  dual  diameter  roller  bearing  concept  for  further  con¬ 
sideration.  This  selection  was  based  upon  lower  cage  speeds  obtained  in 
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the  dual  diameter  bearing  which  result  in  signiiican~ly  iowur  : 
mechanic  losses  (viscous  loss  due  to  shearing  of  the  oil) .  These  lower 
drags  or  losses  in  turn  result  in  a  lower  cage  slip  threshold. 

Table  I 

BEARING  OPERATING  REQUIREMENTS 

Normal  Operating  Temperature 
Operating  Temperature  Range 
Normal  Operating  ON  (mm  x  rpm) 

Operating  Radial  Load  Range  (lb.) 

L | q  Life,  Hr.  (greater  than) 

Radial  Deflection  at  Max.  Load,  (in.) 

Minimum  Heat  Generation 

The  dual  diameter  roller  bearing  configuration  was  subjected  to  a 
thorough  analysis  from  which  a  final  design  was  developed.  A  complete 
systems  analysis  incorporating  the  latest  available  elastohydrodynamic 
lubrication  technology  and  including  fluid  drag  forces  and  thermal  ef¬ 
fects  was  accomplished  by  means  of  a  digital  computer  program  (Reference 
1).  This  analysis  was  used  to  finalize  the  bearing  design  and  determine 
operating  clearances  and  oil  flow  or  lubrication  requirements  of  the 
bearing.  Complete  detailed  working  drawings  of  the  prototype  test  bear¬ 
ings  are  contained  in  the  body  of  this  report. 

Six  prototype  dual  diameter  roller  bearings  were  manufactured  by 
the  Federal-Mogul  Corporation,  Bower  Bearings  Div.,  Detroit,  Michigan, 
to  the  drawings  and  specifications  contained  in  this  report.  A  summary 
of  inspection  data  relating  to  dimensional  quality  of  the  test  bearings 
is  also  included. 

Two  dual  diameter  prototype  roller  bearings  were  tested  by  the 
Midwest  Aero  Industries  Div.,  Pure  Carbon  Co.,  St.  Mary’s,  Pa.  on  a 
special  high  speed  mainshaft  test  rig.  The  test  bearings  were  operated 


600°  F 

-65  to  +800°F 
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initially  with  MIL—L-7808  lubricant  at  250°F  oil  inlet  and  then  later 
with  Folyphenyl  Ether  5P4E  lubricant  at  600°F. 
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SECTION  II 


SELECTION  OF  PRELIMINARY  CONFIGURATION 


1.  CANDIDATE  BEARING  CONFIGURATIONS 

Five  >asic  concepts  or  candidate  bearing  configurations,  as  shown 
on  Figure  l(a-d),  were  investigated.  The  reasons  f>r  considering  the 
various  configurations  are  listed  below: 

a.  Solid  rollers  -  inner  land  riding  cage:  The  land  riding 

cage  solid  roller  concept,  Figure  la,  is  one  of  the  most  com¬ 
monly  used  mainshaft  roller  bearing  designs.  This  design  con¬ 
sists  of  an  inner  race  with  two  integral  shoulders  to  guide 
the  roller  complement.  The  one  piece  fully  machined  cage  has 
close  clearance  and  rides  on  an  oil  film  in  contact  with  the 
lands  on  the  inner  race.  The  outar  race  is  a  sleeve  construc¬ 
tion  which  allows  complete  axial  freedom  of  the  inner  race 

and  roller  complement.  It  is  common  in  some  designs  to  provide 
coolant  flow  through  slots  under  the  inner  race  since  the  inner 
race  operating  temperature  is  approximately  50eF  hotter  than 
the  outer  race  operating  temperature.  An  advantage  of  this 
type  of  cage  construction  is  the  fact  that  friction  in  the  cage 
to  land  areas  will  help  to  drive  the  cage  at  synchronous  or 
epicyclic  speed. 

b.  Solid  roller  -  outer  land  riding  cage:  This  outer  land  rid¬ 
ing  concept.  Figure  lb,  is  very  similar  to  the  concept  above; 
the  main  difference  being  that  the  integral  guiding  flanges 
and  the  cage  to  race  land  guidance  are  on  the  outer  race. 

This  bearing  configuration  is  also  common  in  modern  jet  engines. 
The  advantage  of  this  type  of  construction  is  better  retention 
of  oil  at  the  outer  race  contact  (due  to  centrifugal  forces) 
under  momentary  oil  starvation.  Also  some  engine  applications 
are  more  easily  assembled  with  this  type  of  roller  guidance. 

It  should  be  noted  that  in  this  configuration  the  friction  in 
the  cage  to  land  surfaces  is  a  dissipative  effect  which  will 
tend  to  place  a  drag  on  the  cage  and  add  to  the  cage  slipping 
tendency . 

c.  Hollow  roller  concept:  The  roller  bearing  concept  Is  shown  in 
Figure  lc.  This  concept  can  be  applied  to  either  an  inner 
land  riding  cage  or  an  outer  land  riding  cage  configuration. 
Basically  the  use  of  a  hollow  roller  reduces  the  outer  race  to 
roller  loading.  The  centrifugal  force  on  the  outer  race  due  to 
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ROLLER  DIAMETER 
0.375*“<  D  <  0.625' 


I40MM 
(5.5118") 
BORE  OIAM. 


Figure  1(a).  Solid  Rollers  -  Inner  Land  Riding  Cige 


HOLLOW  ROLLER  CONCEPT 
(Reduce  Roll  Mass) 


the  cage  and  roller  complement  rotative  speed  results  in  signifi¬ 
cant  contact  loading.  Hollow  rollers  of  less  mass  than  solid 
rollers  are  found  to  be  beneficial.  A  note  of  caution  is  indi¬ 
cated  in  that  prior  exoerience  shows  that  too  thin  a  wall  section 
in  the  hollow  rollers  will  lead  to  roller  breakage  due  to  reverse 
bending  of  the  rings.  Thus  mass  reduction  must  be  compromised 
with  ring  bending  fatigue  properties.  Attempts  have  been  made  in 
the  past  to  use  three  or  more  oversized  hollow  rollers  equally 
spaced  in  a  solid  roller  complement  to  provide  rollers  which 
are  always  .nder  load  and  essentially  drives  the  cages  up  to  cyc¬ 
lic  speed.  The  effect  which  was  analyzed  in  the  preliminary 
selection  of  configurations  was  the  mass  reduction  effect 
only. 

d.  Series  hybrid  bearing:  The  series  hybrid  bearing  is  a  relatively 
new  concept.  Figure  Id,  for  reducing  the  effective  speed  of 
rotation  at  the  inner  race  of  a  mainshaft  roller  bearing.  This 
concept  consists  r'  a  rolling  element  bearing  and  an  oil  film 
journal  bearing  in  series.  The  inner  race  of  the  roller  bear¬ 
ing  is  separated  from  the  outer  surface  of  the  rotating  shaft 

by  means  of  the  hydrodynamic  oil  film  bearing.  Both  bearings 
carry  the  full  radial  load,  however,  the  inner  race  of  the 
rolling  bearing  will  rotate  at  some  intermediate  speed  to  the 
shaft  speed  resulting  in  a  lower  effective  DN  operating  regime 
for  the  roller  bearing.  The  sp--ed  reduction  is  the  main  ad¬ 
vantage  of  this  configuration.  Additional  under  race  cooling 
of  the  inner  race  can  also  be  accomplished  if  locating  oil 
film  thrust  bearing  surfaces  are  built  as  pumps  to  supply  ad¬ 
ditional  coolant  flow  under  the  inner  race.  These  locating 
thrust  bearings  to  maintain  the  location  of  the  inner  race  are 
not.  shown  on  Figure  Id,  The  main  disadvantage  is  the  obvious 
one  of  mechanical  complexity. 

e.  Dual  diameter  roller  concept:  The  dual  diameter  roller  con¬ 
cept  is  shown  in  Figure  2.  This  is  basically  a  wagon  wheel 
shaped  roller  with  a  large  central  diameter  and  two  smaller 
diameters  which  look  very  much  like  an  axle  through  the  main 
wheel.  The  purpose  of  this  configuration  is  to  reduce  cage 
operating  speed.  A  reduction  in  cage  orbital  speed  results 

in  a  significant  reduction  of  the  outer  race  to  roller  contact 
load  due  to  centrifugal  effects. 

The  relative  kinematics  of  a  conventional  roller  bearing  and  a  dual 
diameter  roller  bearing  are  shown  in  Figure  3,  The  large  outer  diameter 
of  the  two-diameter  roller  contacts  the  inner  race  surface.  The  email 
roller  diameters  of  the  two-diameter  roller  contact  the  outer  race  sur¬ 
face.  The  effects  of  the  relative  kinematics  are  shown  graphically  in 
Figure  4.  With  a  small  roller  (50%  of  the  major  roller  diameter)  it  can 
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Figure  3,  Kinematics  of  Dual  Dian;,  Roller 
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Figure  4.  Relative  Kinematics  of  Dual  Diameter  Roller 
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be  seen  that  for  the  same  inner  race  speed,  a  cage  rotative  speed  1/3 
less  than  that  of  the  conventional  bearing  can  be  achieved.  At  the  same 
time  the  roller  rotative  speed  about  its  own  axis  is  about  3 /3  greater 
in  the  dual  diameter  roller  than  in  the  conventional  roller  design. 

This  increased  roller  rotative  speed  enhances  the  entrainment  velocity 
which  is  important  in  the  development  of  a  full  elastohydrodynamic  oil 
film  in  the  contacts.  Higher  entrainment  velocities  of  the  two-diameter 
roller  will  essentially  result  in  thicker  EHD  oil  films  in  the  contact. 
Thus  a  4.0  x  10°  DN  bearing  of  the  two-diameter  roller  configuration  will 
result  in  cage  rotative  speeds  and  resulting  centrifugal  effects  at  the 
roller  to  outer  race  contact  comparable  to  those  found  in  a  2.67  x  10^ 

DN  bearing  of  conventional  design. 

2.  LOAD-LIFE  PARAMETRIC  STUDY 

Initial  calculations  indicated  that  the  requirements  for  greater 
than  1,000  hours  L-10  operating  life  at  4.0  x  10^  DN  could  not  be 
achieved  with  2,000  lbs,  radial  load.  The  2,000  lb,  load  is  repre¬ 
sentative  of  momentary  maneuvering  loads  and  not  intended  as  a  long 
time  operating  load  condition  for  design  of  mainshaft  bearings.  The 
duty  cycle  given  in  Table  II  was  used  for  the  parametric  study. 

A  special  short  length  computer  program  was  used  for  the  parametric 
study  and  consisted  of  a  full  load  deflection  analysis  of  the  bearing 
with  elastohydrodynamic  films  in  the  race  roller  contacts.  The  EHD 
film  thicknesses  were  computed  per  the  classical  isothermal  Dowson- 
Higginson  formulas  (Reference  2) .  This  is  a  simplification  in  that  no 
cage  slip  was  considered »  However,  for  the  purposes  of  load  life  studies 
this  is  a  valid  assumption  especially  for  preliminary  configuration 
selection.  The  three  radial  loads  identified  in  Table  II  were  used  to 
analyze  the  operating  L-10  life  and  the  results  were  prorated  over  the 
percentages  of  operating  times  (Reference  3) .  The  figures  presented 
are  in  terms  of  the  prorated  life  over  the  duty  cycle , 
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Table  II 


DUTY  CYCLE  AT  4.0  x  106  DN 


%  Time 

Radial  Load 

10 

2000  pousiJs 

20 

1000  pounds 

70 

500  pounds 

The  parametric  studies  were  run  with  the  following  assumptions: 

•  Constant  pitch  diameter  =6.50  inches  =  E 

•  Number  of  rollers  =  tt  times  E  divided  by  1.5  times  D 

•  Effective  length  cf  contact  =  roller  diameter  (D) 

«  4  x  10^  DN  (28,740  shaft  rpra  at  140  mm  bore) 

a  Dual  diameter  roller  .8D  outer  length,  ,5D  inner  length 

The  L-10  prorated  bearing  life  over  the  entire  duty  cycle  in  hours 
for  the  hollow  roller  concept  is  shown  in  Figure  5  where  bearing  life  is 
plotted  as  a  function  of  roller  diameter.  The  bottom  curve  on  Figure  5 
represents  a  solid  roller  of  the  configurations  described  in  Figure  la 
and  lb.  A  similar  plot  of  the  L-10  prorated  bearing  life  over  the  duty 
cycle  for  the  dual  diameter  roller  bearing  is  shown  in  Figure  6.  An 
obvious  knee  in  the  curves  is  apparent  for  all  of  the  roller  concepts 
at  a  roller  diameter  of  0.5  inches.  This  can  be  explained  in  terms  cf 
the  exponential  reduction  in  life  with  increasing  load.  The  larger 
diameter  rollers  have  a  higher  outer  race  contact  load  because  of  in¬ 
creased  centrifugal  force  and  over  0.5  inch  roller  diameter  this  increase 
in  centrifugal  effects  is  greater  than  tne  corresponding  increase  in 
dynamic  capacity  of  the  roller  race  contact. 

The  results  of  the  parametric  study  in  terms  of  the  design  goal  of 
more  than  1,000  hours  L-10  life  are  (a)  roller  diameter  must  be  0.5  inches 
or  less  for  all  configurations,  (b)  solid  rollers  are  net  acceptable  in 
any  diameter  size,  (e)  hollow  rollers  should  be  0.5  to  0.7  hollow,  (a)  dual 
diameter  rollers  with  0.5  to  0.6  small  to  large  diameter  roller  ratio 
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are  of  interest,  (e)  the  requirement  for  0.01  inches  maximum  radial  de¬ 
formation  under  2,000  lbs.  load  indicates  that  hollow  rollers  over  0.7 
hollowness  are  not  recommended.  The  parametric  life  study  also  shows 
that  L-.IO  fatigue  life  can  tve  a  limiting  factor  at  4.0  x  10^  DN 
operation. 

3.  CAGE-SLIP  AND  HORSEPOWER  LOSS  STUDY 

The  most  promising  configurations  resulting  from  the  parametric 
study  were  used  in  a  full  EHD  computer  analysis  (Reference  1)  containing 
the  latest  available  elastohydrodynamic  (EHD)  technology  (Reference  4). 

The  majority  of  the  computer  runs  were  made  with  MIL-L-7808  oil  at  250°F 
oil  inlet  temperature.  One  check  run  each  for  a  single  diamet*  and 
a  dual  diameter  roller  configuration  were  made  with  Polyphenyl  Ether 
5P4E  lubricant  at  500°F  oil  inlet  temperature.  The  actual  run  data  is 
summarized  in  Table  III,  which  also  contains  a  summary  of  results  includ¬ 
ing  percentage  cage  slip  and  horsepower  dissipation.  More  detailed 
information  is  summarized  in  Table  XV.  L-10  life,  roller  load,  film 
thickness,  and  specific  film  thickness  data  as  well  as  sliding  speeds 
at  both  the  outer  and  inner  races  are  given  for  the  maximum  loaded  roller. 

The  primary  problem  encountered  in  high  speed  malnshaft  roller 
bearings  is  cage  slip.  This  is  the  phenomena  which  results  in  high 
sliding  velocities  between  the  rollers  and  the  inner  race  contacts. 

Severe  glazing  and  micro-pitting  results  in  bearing  failure.  The  full 
EHD  computer  analysis  gives  detailed  information  regarding  slip  condi¬ 
tions.  Figure  7  shows  the  percentage  cage  slip  as  a  function  of  radial 
load  for  various  operating  conditions.  An  inner  race  land  riding  cage 
with  hollow  rollers  has  considerable  slip  and  is  used  as  a  reference 
point  in  Figure  7.  Seven  hundred  pounds  radial  load  are  required  for 
non-slip  cr  synchronous  running  of  the  cage.  The  outer  race  land  riding 
hollow  roller  design  has  considerably  higher  slip  than  any  of  the  other 
configurations  and  only  one  check  run  is  plotted  on  Figure  7.  The  bene¬ 
ficial  effects  of  having  the  frictional  torque  and  the  cage  land  con¬ 
tact  acting  as  a  driving  torque  indicates  that  an  inner  land  riding  cage 
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SCHEDULE  OF  COMPUTER  RUNS 
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configuration  should  be  considered.  For  this  reason  all  dual  diameter 
roller  computer  runs  were  made  with  an  inner  land  riding  cage  configura¬ 
tion.  The  dual  diameter  roller  bearing  has  significantly  less  cage 
slip  as  shown  on  Figure  7.  Another  way  to  represent  sliding  or  cage  slip 
is  to  show  the  reduction  in  cage  rpm  as  a  function  of  radial  load. 

This  is  given  in  Figure  8.  It  is  immediately  apparent  that  the  hollow 
roller  inner-land-riding  cage  speed  or  a  single  diameter  roller  bearing 
concept  cage  speed  is  considerably  higher  than  the  corresponding  cage 
speed  for  the  dual  diameter  roller  concept.  As  evident  in  Figure  8,  a 
much  lower  load  is  required  to  prevent  sliding  or  iage  slip  in  the  dual 
diameter  roller  concept. 

An  object  of  the  present  study  is  to  minimize  heat  generation  in 
the  bearing.  An  indication  of  heat  generation  is  the  computed  horse¬ 
power  loss  of  the  bearing.  As  shown  on  Figure  9,  the  horsepower  losses 
for  both  the  inner  land  riding  cage  configurations  (single  diameter 
roller  and  dual  diameter  roller)  are  essentially  the  same.  Horsepower 
losses  are  expected  to  be  somewhat  lower  with  the  use  of  Polyphenyl 
Ether  5P4E  lubricant  as  compared  to  MIL-L-7808  oil.  The  oil  viscosity 
of  the  5P4E  at  500°F  is  less  than  that  of  the  7808  at  250°F.  Viscous 
drag  losses  predominate  in  the  bearing  at  4.0  x  10  DN.  Therefore,  lower 
oil  viscosity  (at  synchronous  cage  speed)  will  result  in  less  power  loss. 

4.  SERIES  HYBRID  BEARING 

The  series  hybrid  bearing  concept,  Figure  Id,  consists  of  an  oil 
film  bearing  between  the  inner  race  of  the  roller  bearing  and  the  shaft. 
The  concept  is  one  of  speed  sharing  or  speed  reduction.  The  rotative 
speed  of  the  roller  bearing  inner  race  is  determined  by  the  frictional 
torque  of  the  oil  film  journal  bearing  matching  the  frictional  torque 
of  the  roller  bearing.  The  advantage  of  speed  sharing  is  the  fact  that 
reduction  speed  on  the  inner  race  of  the  roller  bearing  results  in  a 
lower  operating  PN  speed  regime. 

A  three  lobe  oil  film  hydrodynamic  bearing  was  considered  for  this 
application.  Satisfactory  bearing  stability  (absence  of  half-frequency 
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Figure  9.  Horsepower  Loss  vs.  Speed 


whirl)  can  be  obtained  with  this  type  of  bearing.  Good  inherent 
stability  was  shown  to  exist  under  a  shaft  weight  of  35  lbs.  (represent¬ 
ing  an  anticipated  test  rig  shaft).  This  stability  decreases  as  shaft 
weight  increases,  however,  satisfactory  stability  for  the  bearing  under 
study  could  be  obtained  up  to  200  lbs.  shaft  weight.  Minimum  oil  films 
on  the  order  0.001  inch  can  be  anticipated  under  light  200  lb.  loadings. 
The  oil  films  would  reduce  to  0.0004  inches  under  the  maximum  2,000 
lbs.  applied  external  load.  The  frictional  torque  of  a  three-lobe  jour¬ 
nal  bearing  was  computed.  The  roller  bearing  inner  race,  Figure  Id, 
when  supported  on  an  oil  film  journal  bearing  only,  would  be  completely 
free  to  move  axially  with  respect  to  the  shaft.  Thus  locating  oil  film 
bearings  would  be  required  to  maintain  the  position  of  the  roller  bear¬ 
ing  inner  race.  Locating  thrust  loads  would  be  small.  One  approach 
is  to  use  a  conical  type  bearing  (Reference  5) .  A  pair  of  thrust  bear¬ 
ings  operating  against  the  face  of  the  roller  bearing  inner  race  were 
envisioned  as  a  series  of  step  pad  thrust  surfaces.  The  frictional 
torque  of  such  surfaces  amounted  to  approximately  15  or  20%  of  the  over¬ 
all  journal  bearing  friction.  This  amount  was  added  to  the  oil  film 
bearing  system  torques. 

The  operation  of  the  journal  bearing  would  be  in  the  turbulent 
regime  and  factors  for  the  computation  of  turbulent  friction  are  usu¬ 
ally  applied  to  laminar  friction  estimates.  This  approach  was  used  in 
arriving  at  the  total  journal  bearing  friction  curve.  The  friction  loss 
vs.  speed  ratio  for  MIL-L-78C8  oil  at  250°  is  shown  in  Figure  10.  The 
solid  lines  indicate  the  friction  torque  expressed  in  inch  pounds  of  the 
oil  film  journal  and  thrust  bearing  surfaces  is  a  function  of  the  ratio 
of  the  roller  bearing  inner  race  speed  to  shaft  speed.  The  dashed  curve 
represents  the  frictional  torque  of  the  inner  land  riding  cage  hollow 
roller  bearing  configuration.  The  point  of  intersection  of  these  two 
frictional  torque  curves  defines  the  operating  speed  of  the  roller  bear¬ 
ing  inner  race.  It  is  seen  in  Figure  10  that  the  present  bearing  system 
would  result  in  a  speed  ratio  of  0.84.  Recent  work  by  NASA  (Reference 
5)  with  a  series  hybrid  concept  using  a  ball  thrust  bearing  and  an  oil 


Figure  10.  Friction  Loss 
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film  thrust  bearing  surface  resulted  in  a  speed  sharing  ratio  of  ap¬ 
proximately  0.66  at  best.  The  analysis  is  highly  dependent  upon  accurate 
calculation  of  frictional  torques.  The  frictional  torques  even  in  the 
turbulent  regime  of  the  fluid  bearing  has  been  fairly  well  established 
and  correlated  with  experimental  results  and  can  be  considered  realistic. 
High  speed  roller  bearing  torque  needs  to  be  carefully  verified  to  de¬ 
termine  the  anticipated  speed  sharing  ratio.  It  is  possible  to  add 
frictional  torque  to  the  roller  bearing  concepts  by  increasing  the  amount 
of  oil  forced  through  the  bearing.  Thus,  it  is  possible  by  raising  the 
frictional  torque  of  the  roller  bearing  to  obtain  a  more  favorable  speed 
sharing  ratio  as  could  be  seen  by  raising  the  curve  on  Figure  10.  How¬ 
ever,  any  reductions  in  roller  bearing  torque  would  result  in  reduction 
in  the  anticipated  speed  sharing  ratio. 

5.  SELECTION  OF  DUAL  DIAMETER  ROLLER  BEARING 

The  dual  diameter  roller  bearing  concept,  Figure  2,  with  an  inner 
land  riding  cage  was  recommended  for  the  thorough  design  analysis, 
fabrication  and  prototype  testing.  Selection  of  the  dual  diameter  rol¬ 
ler  bearing  concept  was  based  upon  the  lower  cage  speeds  obtained  which 
result  in  significantly  less  cage  slip. 

The  series  hybrid  bearing  concept  is  interesting  but  was  not  recom¬ 
mended  for  further  study  under  the  present  program.  Operation  of  the 
series  hybrid  bearing  design  is  highly  dependent  upon  accurate  determina¬ 
tion  of  roller  bearing  friction  torque.  The  horsepower  losses  of  both 
conventional  and  dual  diameter  roller  bearing  concepts,  Figure  9,  are 
3hown  to  be  essentially  the  same.  Therefore  the  selection  of  a  roller 
bearing  configuration  for  further  development  at  this  time  does  not 
affect  future  consideration  of  development  of  the  series  hybrid  bearing 
concept. 
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SECTION  III 


LUBRICANT  SELECTION 

A  summary  of  full  scale  bearing  test  results  of  8  candidate  high 
temperature  (600“F)  lubricants  is  given  in  Table  V,  which  was  taken  from 
Reference  6.  All  of  the  lubricants  in  this  study  were  operated  with  an 
inerting  blanket  of  nitrogen  gas.  Another  recent  study  (References  7 
and  8)  considered  five  candidate  high  temperature  lubricants.  In  this 
case  Table  VI,  all  of  the  lubricants  were  again  protected  by  an  inert 
nitrogen  blanketing  ga9  with  the  exception  of  Polyphenyl  Ether  5P4E. 

The  Polyphenyl  Ether  was  operated  in  an  oxygen  atmosphere.  The  only 
advanced  lubricant  showing  reasonable  success  without  special  inert 
oxidation  protection  is  the  5P4E  Polyphenyl  Ether.  Figure  11  is  a 
reproduction  of  a  report  of  bearing  fatigue  life  experiments  with  Poly¬ 
phenyl  Ether  (Reference  8) .  Indications  are  that  Pclyphenyl  Ether  is 
an  acceptable  high  temperature  lubricant  up  to  600° F  but  that  the  failure 
mode  may  be  one  of  moderate  glazing  and  micropit ting  of  the  races  and 
rolling  elements.  This  is  attributed  to  thinner  EH")  lubricant  films  in 
the  operating  race  and  roller  contacts.  The  various  bearing  configura¬ 
tions  were  analyzed  with  5P4E  Polyphenyl  Ether  at  500°F  to  determine  the 
operating  specific  film  thicknesses  in  the  roller  race  contacts. 

The  properties  of  MIL-L-7808  and  Polyphenyl  Ether  lubricants  are 
discussed  in  the  elastohydtodynamic  lubrication  preliminary  design 
manual  (Reference  4) .  Additional  data  (Reference  9)  regarding  high  pres¬ 
sure  viscometer  experiments  witn  Polyphenyl  Ether  5P4E  were  also  used 
for  estimating  properties  of  this  lubricant  at  500°F.  Lubricant  data 
as  a  function  of  operating  temper?  ure  is  given  in  Table  VII  for  the 
MIL-L-7808  aud  5P4E  lubricants.  This  ’ita  as  a  function  of  temperature 
is  consistent  with  the  extrapolation  routines  of  Reference  4. 


SUMMARY  OF  TEST  RESULTS 
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BEARING  lift,  MILLIONS  OF  INNER -RACE  REVOLUTIONS 


Figure  11.  Bearing  Life  with  5P4E  Lubricant 

Rolling-element  fatigue  life  of  120-mm  bore  angular- 
contact  ball  bearings  run  with  polypheny!  ether  lubricant. 
Material,  AISI  M-50  steel;  thrust  load,  4365  lb;  speed, 
12,000  rpm;  temperature,  600°F  air  environment;  [failure 
index,  2  out  of  26.] 
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SECTIOfl  IV 


DUAL  DIAMETER  SEARING  ANALYSIS  AND  DESIGN 
1.  DESIGN  CONSIDERATIONS 

The  dual  diameter  roller  bearing  concept  was  selected  as  a  result 
of  the  preliminary  configuration  study.  Selection  was  based  upon  the 
lower  cage  speeds  obtained  which  will  result  In  significantly  less  cage 
slip.  The  orbital  cage  speed  of  the  dual  diameter  roller  bearing  with 
a  diameter  ratio  of  0.5  will  result  in  a  1/3  reduction  of  cage  speed  as 
seen  by  Figure  A.  A  major  roll  diameter  of  0.5  inches  was  selected  as 
an  optimum  basid  on  Figure  6.  The  resulting  estimated  bearing  life  (under 
the  duty  cycle  ot  Table  II)  would  be  5,000  hours  L-10  life  shown  on  Figure  6. 
This  is  well  in  excess  of  the  design  goal  of  1,000  hours. 

The  primary  problem  encountered  in  high  speed  mainshaft  roller  bear¬ 
ings  is  cage  slip.  This  is  a  phenomena  which  results  in  high  sliding 
velocities  between  rollers  and  inner  race  contacts.  Severe  glazing  and 
micropitting  results  in  bearing  failure.  The  minimum  radial  load  to 
prevent  skid  or  cage  slip  in  the  dual  diameter  roller  bearing  is  333  lbs. 
as  shown  on  Figures  7  and  8.  This  minimum  radial  load  was  used  in  the 
design  and  performance  evaluation  calculations.  Heavier  radial  loads  will 
result  in  reduced  bearing  fatigue  life,  but  do  not  endanger  tne  bearing 
from  a  skid  or  cage  slip  standpoint  and  are  of  relatively  minor  interest. 

It  is  essential  to  assure  that  cage  slip  will  not  occur  under  operating 
conditions.  A  complete  systems  analysis  (Reference  1)  incorporating  the 
latest  available  elastohydrodynamic  lubrication  technology  (Reference  A) 
and  including  fluid  drag  forces  and  thermal  effects  were  selected  as  the 
best  approach  to  achieve  a  practical  design. 
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2.  INNER  RACE  STRESSES 


a  study  of  the  effect  of  inner  race  thickness  and  total  interference 
fit  between  the  inner  race  and  shaft  (with  the  hollow  test  shaft)  to  deter¬ 
mine  race  circumferential,  hoop  or  tangential  tension  stresses  was  conducted 
using  standard  ring  formulas.  Race  radial  thickness  over  0.25  inches  has 
little  effect  upon  reducing  circumferential  tangential  stresses.  Race 
thicknesses  less  than  0.250  inches  radial  thickness  can  result  in  unacceptably 
high  stresses.  Therefore  the  inner  race  thickness  should  meet  or  exceed 
0.25  inches. 

High  shaft  speeds  (up  to  28,600  rpm)  will  create  significant  circumfer¬ 
ential,  hoop  or  tangential  stresses  in  the  inner  race.  Figure  12  shows  a 
magnitude  of  this  maximum  tangential  stress  as  a  function  of  shaft  speed 
expressed  both  in  rpm  and  DN  values  for  the  140  mm  bore  bearing.  Approx¬ 
imately  50,000  psi  tangential  stress  can  be  anticipated  at  4.0  x  10^  DN 
with  the  0.25  inch  thick  race  of  the  prototype  bearing  design. 

Approximate  physical  properties  of  M-50  fully  hardened  steel  are  given 
in  Table  VIII.  The  material  is  brittle  at  room  temperature  and  care  must  be 
exercised  to  prevent  over-stressing.  The  material  does  become  more  ductile 
as  temperature  increases.  The  effect  of  stresses  due  to  speed.  Figure  12, 
and  stresses  due  to  interference  fit  are  essentially  additive.  Approximately 
0.012  inches  diametral  Interference  fit  would  be  required  at  room  temperature 
to  assure  a  working  interference  fit  of  0.001  at  4.0  x  10  DN  operation. 

This  heavy  interference  fit  would  result  in  almost  60,000  psi  tangential 
ring  stresses  initially.  It  was  decided  to  compromise  with  a  lesser  inter¬ 
ference  fit  at  room  temperature  (0.0063  inches)  and  allow  the  inner  race  to 
become  loose  with  respect  to  the  shaft  at  approximately  3.0  x  10^  DN.  Key 
slots  in  the  inner  race  allow  the  race  to  be  axially  clamped  and  prevent 
relative  rotation  of  the  race  and  shaft.  The  loose  fit  between  "nner  rac_ 
and  shaft  at  speeds  over  3.0  x  10  DN  will  keep  the  Inner  race  tangential 
stresses  to  less  than  60,000  psi. 

High  tensile  ring  stresses  may  have  a  detrimental  effect  on  inner  race 
rolling  contact  fatigue.  Thermal  gradients  may  also  cause  the  interference 


0  7,200  14,300  21,400  28,600  RPM 

SHAFT  SPEED 

Figure  12.  Inner  Race  Maximum  Tangential  Stress  vs.  Shaft  Speed 
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TABLE  VIII 


PROPERTIES  OF  M-50  HARDENED  TO  RC-60 
Per  VASCO  Technical  Bulletin 


Property 

Room  1  snip. 

800°F 

Tensi le  ULT.  (psl) 

411,500 

345,000 

Tens! le  Yield  (psi ) 

388,000 

262,500 

Elongation  (1.5"  gauge)  (%) 

2 

6 

Red.  of  Area  {%) 

2 

17 

35 


fit  to  increase  (depending  upon  the  amount  of  under  inner  race  cooling). 
Therefore  every  effort  was  made  to  limit  inner  race  tangential  stresses 
at  maxim  _u  and  500°F  operation.  A  partially  hollow  or  compliant 

inner  race  sign  Figure  13  was  selected  to  provide  access  for  the  under 
and  through  inner  race  lubrication  scheme.  The  beam  bending  compliancy 
may  also  provide  relief  from  thermal  lockup  or  thermal  preloading  at  high 
temperature  and  high  speed  operation. 

An  existing  FIRL  finite  element  computer  program  was  used  to  evaluate 
the  tangential  ring  stresses  in  the  actual  inner  race  design  at  4.0  x  10^ 
and  3.5  x  10^  BN.  Figure  13  illustrates  the  inner  race  model  and  the 
location  of  various  tangential  stresses.  The  stress  values  of  Figure  33 
are  in  good  agreement  with  the  approximate  solutions  used  to  obtain  the 
values  in  Figure  12.  It  is  readily  apparent  that  stresses  on  the  order 
of  57,000  psi  will  occur  in  the  inner  race  at  4.0  x  10^  DN.  Stress  level 
will  drop  to  approximately  44,000  psi  maximum  at  3.5  x  10°  DN  (25,000  rpm) . 

A  third  case  was  analyzed  where  circumferential  line  loads  were  applied  in 
the  vicinity  of  the  roller  inner  race  contacts  to  simulate  stress  conditions 
under  500  lb.  external  radial  load.  It  can  be  seen  from  Table  IX  that  the 
stress  values  range  between  a  high  of  43,000  psi  and  a  low  of  29,000  psi. 

The  only  available  data  to  date  for  M-50  steel  at  500°F  in  reverse  bending 
f,  :igue  (Reference  10)  indicates  a  life  which  is  asymptotic  at  60,000  psi 

g 

b'-'ding  stresses  beyond  10  cycles.  The  fluctuation  of  inner  race  bending 
stresses  are  all  unidirectional  and  remain  tensile  as  shown  in  Table  IX 
and  do  not  exceed  a  maximum  value  of  43,000  psi.  The  inner  race  compliant 
section  calculates  to  have  adequate  reverse  bending  fatigue  properties. 

3.  LUBRICATION 

The  method  of  lubrication  for  the  dual  diameter  roller  bearing  is 
Illustrated  in  Figure  14.  Jet  oil  is  introduced  to  an  undercut  in  the 
test  shaft  as  shown.  Oil  then  passes  through  the  shaft  to  the  undercut 
portion  of  the  inner  race  whicl  'cts  as  an  oil  reservoir.  Centrifugal 
force  resulting  from  shaft  rotation  then  provides  pressure  to  force  a 
flow  of  oil  through  orifice  holes  in  the  inner  race.  This  oil  is  then 
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Figure  13.  Finite  Element  Analysis  of  Inner  Race 
Tangential  Stresses  at  4.0  x  10&  DN 


TABLE  IX 


INNER  RACE  TANGENTIAL  STRESS  (KSI) 


Loading/Element  No. 

1 

10 

81 

90 

171 

208  215 

270 

370 

402 

392 

4.0  x  106  DN 

55 

55 

55 

53 

55 

57 

52 

57 

51 

55 

51 

3.5  x  106  DN 

42 

41 

42 

41 

41 

43 

40 

44 

39 

43 

39 

3.5  x  106  DN 

500  lbs. 

37 

23 

36 

33 

37 

40 

36 

41 

38 

43 

39 

Radial  Load 
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Figure  14.  Thermal  Nodes  ar.i  Methods 


introduced  to  the  underside  of  the  cage  and  is  further  circulated  through¬ 
out  the  bearing  by  the  pumping  action  of  the  rollers.  Auxiliary  jet  oil 
lubrication  into  the  sides  or  faces  of  the  bearing  can  also  be  provided 
as  necessary  to  obtain  satisfactory  cooling. 

Three  computer  analysis  runs  were  executed  using  three  different  flow 
rates  (6.1,  12.0  and  18.0  lbs. /min.)  of  5P4E  at  500°F  oil  inlet  temperature. 
The  minimum  load  to  prevent  skid  of  333  lbs.  and  the  top  design  shaft  speed 
of  28,560  rpm  (4.0  x  10^  DN)  were  used.  A  complete  thermal  analysis  was 
included  in  the  solutions.  A  summary  of  the  effects  of  oil  flow  upon 
bearing  performance  is  given  in  Table  X  and  shown  graphically  in  Figure  15. 
Examination  of  Figure  15  shows  no  clear  maxima  or  minima  characteristics. 

It  is  evident  that  increased  oil  flow  rate  results  in  reduced  outlet  tempera¬ 
ture  and  increased  bearing  horsepower  loss.  The  outlet  oil  temperature 
shows  some  indication  of  leveling  out  below  6  lbs.  per  minute,  but  at 
unacceptably  high  values  in  excess  of  700°F.  Increasing  the  oil  flow 
from  6  to  12  lbs.  per  minute  drops  the  outlet  oil  temperature  by  51.8°F. 
Increasing  the  oil  flow  from  12-18  lbs.  per  minute  drops  the  outlet  oil 
temperature  by  only  22.3°F. 

An  oil  flow  rate  of  12  lbs. /min.  was  selected  as  the  design  point. 

This  flow  can  be  obtained  with  6  orifices  of  0.028  inches  diameter.  The 
calculated  oil  outlet  temperature  of  4.0  x  10^  DN  is  598. 5°F  which  is 
consistent  with  tht  program  goal  of  600°F  operation.  Table  XI  gives  the 
resulting  oil  flow  rates  for  5P4E  at  50G°F  oil  inlet  and  MIL-L-7808  at 
250°F  oil  inlet  as  a  function  of  shaft  speed. 

4.  PERFORMANCE  ANALYSIS 

Eight  computer  solutions  were  executed  at  various  shaft  speeds  with 
5P4E  and  MIL-L-7808  oil  lubricant  under  a  constant  bearing  radial  load 
of  333  pounds.  A  summary  of  these  computer  runs  and  performance  results 
are  contained  in  Table  XII.  Initial  fit-up  conditions  are  0.0063  inches 
shrink  fit  between  bearing  inner  race  and  the  shaft  and  0.011  inches  initial 
(machined  or  ground)  diametral  clearance.  All  of  the  solutions  summarized 
in  Table  XII  included  a  full  thermal  analysis  consisting  of  21  thermal  modes 
described  in  Figure  14.  Description  of  the  node  locations  and  the  temperature 
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TABLE  X 

EFFECT  OF  OIL  FLOW  RATE  ON  BEARING  PERFORMANCE 


Run  2P 

3P 

4p 

Flow  (lbs/min) 

6.1 

12.0 

18.0 

HP  (loss) 

21.7 

23.53 

25.4 

Temp.  Out 
of  Brg.  (°F) 

650.3 

598.5 

576.2 

Sump.  Temp.  (*F) 

608.1 

561.8 

544.4 

Shaft  Speed  28,560  RPM  (4.0 
Lubricant;  5P4E  @500°F  oil 
Radial  Load:  333  lbs 

x  10&  DN) 
inlet 

TABLE  XI 

OIL 

FLOW  RATE  VS. 

SPEED  (DN) 

Lubricant 

Inlet  Temp.  °F 

4.0xl06  DN 

Oil  Flow  Rate  Ibs/min 
3.0xl06  DN  2.0xl06  DN 

l.OxlO6  DN 

Polypheny! 
Ether  5P4E 

500 

12.0 

9.03  6.05 

3.0 

MIL-L-7808 

250 

8.98 

6.82  4.60 

2.33 

Six  0.028" 

Diam.  Orificies  in  Inner  Race 
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Figure  15.  Effect  of  Oil  Flow  Rate  Upon  Outlet 
Temperature  and  Bearing  Power  Loss 


profiles  for  the  8  performance  runs  are  summarized  in  Table  XIII. 

The  temperatures  when  using  the  Polyphenyl  Ether  5P4E  lubricant  are 
acceptable  at  all  speeds.  The  average  oil  temperature  within  the  bearing 
(node  11)  at  4.0  x  10^  DN  operation  of  631.9°F  is  acceptable  for  short 
duration  test  runs.  Operating  temperatures  with  MIL-L-7808  oil  are  not 
acceptable  and  auxiliary  side  jet  lubrication  was  required  during  the 
initial  checkout  test  runs  with  this  lubricant. 

The  predicted  frictional  horsepower  losses  are  shewn  in  Figure  16  as 
a  function  of  the  speed  parameter  DN.  Horsepower  losses  are  also  summarized 
in  Table  XII.  A  total  loss  of  22.6  horsepower  was  predicted  at  4.0  x  10  DN 

using  Polyphenyl  Ether  5P4E  at  5006F  oil  inlet  temperatures. 

Radial  displacements  are  summarized  in  Table  XII  and  include  one  half 
of  the  operating  internal  clearance  in  the  bearing.  The  specifications 
require  that  the  total  radial  displacement  not  exceed  0.010  inches  under 
2000  lbs.  radial  load.  An  additional  elastic  deflection  between  rollers 
and  raceways  of  0.0006  inches  should  be  added  to  the  values  in  Table  XII 
to  adjust  for  the  difference  between  the  333  lb.  analysis  load  and  the 
2000  lb.  specific  tion  load.  All  total  displacement  values  are  well  below 
(less  than  half)  maximum  allowable  values.  The  largest  radial  displacements 

occur  at  the  lowest  shaft  speeds.  The  effects  of  shaft  speed  and  thermal 

effects  in  the  bearing  tend  to  decrease  diametral  clearance.  A  summary  of 
the  fluid  film  forces  and  roller  contact  forces  acting  on  the  cage  are 
given  in  Table  XIV.  It  is  evident  from  the  table  that  the  drwing  torque 
resulting  from  the  loaded  rollers  contacting  the  cage  pockets  is  not 
sufficient  to  prevent  serious  cage  slip  at  4.0  x  10^  DN  operation.  The 
fluid  film  drag  on  the  inner  surface  and  inner  race  to  cage  land  is  necessary 
to  drive  the  cage  at  speed  and  prevent  cage  slip.  An  outer  riding  cage 
with  resultant  fluid  drag  torques  at  the  outer  race  land  contacts  would 

g 

result  in  significant  cage  slip  at  4.0  x  10  DN.  Roller  to  cage  web  contact 
forces  of  10.3  pounds  are  predicted  at  the  maximum  DN  operation. 
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5PAE  Input  Temperature  500° F 
7808  Input  Temperature  250° F 
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(2)  Has  5  rollers  in  contact 
All  others  have  three. 


5.  EHD  ^METERS 


The  cter  race  to  roller  contact  EHD  conditions  for  unloaded  rollers 
are  sumrnai  ed  in  Tabic  XV.  The  outer  race  to  roller  conditions  for  the 
maximum  loaded  roller  in  each  case  are  summarized  in  Table  XVI.  The  inner 
race  contact  EHD  conditions  for  the  maximum  loaded  roller  are  summarized 
in  Table  XVII. 

Thin  elastohydrodynamic  (EHD)  oil  films  (less  than  3.0  x  10  ^  inches) 
are  predicted  with  Polyphenyl  Ether  5P4E  lubricant  at  500°F  oil  inlet  tempera¬ 
ture  for  the  bearing.  The  calculated  oil  inlet  temperature  at  the  entrance 
to  the  EHD  contacts  are  given  in  the  appropriate  tables.  The  thin  EHD 
films  result  in  low  (approximately  0.25)  specific  film  values.  Specific 
film  thickness  is  the  ratio  of  the  actual  oil  film  thickness  to  the  root 
mean  square  of  the  contacting  surface  finishes.  Operation  at  600°F  or 
higher  temperatures  with  5P4E  lubricant  can  be  expected  to  result  in  surface 
oriented  damage  or  distress  prior  to  any  evidence  of  classical  sub-surface 
oriented  fatigue  damage. 

Somewhat  thicker  EHD  films  are  predicted  with  MIL-L-7808  oil  applied 
to  the  bearing  at  250°F.  The  specific  films  are  approximately  1.0  which 
is  still  in  the  surface  related  distress  region  of  operation. 

The  sliding  contribution  to  the  total  traction  values  between  roller 
and  race  contacts  is  predominant  for  both  contacts  of  the  maximum  loaded 
rollers  (Tables  XVI  and  XVII).  The  rolling  contribution  to  the  total 
traction  is  significant  at  the  outer  race  contact  of  the  unloaded  rollers. 
Table  XV.  All  of  the  rolling  contacts  are  in  the  low  slip  (sliding  velocities 
less  than  1.0  inches  per  second)  region,  with  the  exception  of  1.0  x  10  DN 
operation  with  5P4E  lubricant. 

6.  DESIGN  DRAWINGS 

The  five  drawings  (Figures  17-21)  completely  describe  the  final  design 
of  the  dual  diameter  roller  bearing. 
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.  Two  Diameter  Roller  Bearing  Assembly 


Diameter  Roller  Bearing 
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Figure  20.  Two  Diameter  Roller 
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SECTION  V 


DUAL  DIAMETER  BEARING  FABRICATION 

Six  prototype  test  bearings  Figure  22  were  fabricated  to  the  design 
drawings  (Figures  17-21)  by  the  Bower  Bearing  Div.,  Federal-Mogul  Corporation, 
Detroit,  Mich.  A  summary  of  the  six  prototype  bearing  .dimensional  sizes 
and  variations  are  given  in  Table  XVIII. 

The  inner  race  roller  path  circularity  and  contour  for  test  bearing 
No.  3  are  shown  in  Figures  23  and  24  and  are  representative  of  all  of  the 
bearings.  The  circularity  and  raceway  contour  of  the  outer  race  roller 
path  are  shown  in  Figures  25  and  26  for  te3t  bearing  No.  3. 

The  roller  design  for  the  dual  diameter  roller.  Figure  20,  is  not 
an  easy  shape  to  produce.  However,  the  test  bearing  rollers  were  centerless 
ground  and  crowns  were  maintained  within  drawing  limits  on  all  three  roller 
portions,  as  shown  by  the  typical  roller  traces.  Figure  27. 

All  material  specifications,  silver  plating,  balance,  etc.,  are  given 
in  the  design  drawings.  Figures  17-21.  All  prototype  bearings  met  or  exceeded 
the  design  drawing  requirements. 
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Figure  23.  Bearing  Ho.  C3  Inrier  Race  Circularity 


Chant  tiC.  i3;«4 


i igure  24.  Bearing  Ho.  C3  Inner  Raceway  Contour 


ure  26.  Bearing  No.  C3  Outer  Raceway  Contour 


Figure  27.  Typical  Dual  Diameter  Roller 
Crown  Profile  Traces 


SECTION  VI 


DUAL  DIAMETER  BEARING  TESTS 


1 

1 

1 
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1.  TESY  RIG 

The  dual  diameter  bearing  testing  was  performed  at  Midwest  Aero 
Industries,  Div.  Pure  Carbon  Co.,  St.  Mary’s,  Pa.,  on  an  existing  250  hp 
drive  test  stand  with  an  eddy  current  clutch  and  speed  increaser  gear  box. 

The  test  rig  is  shown  in  Figures  28  and  29.  A  cross-sectional  drawing  of 
the  rig  is  given  in  Figure  30.  The  numbered  locations  correspond  to  the 
thermocouple  temperature  readings  of  the  actual  test  data.  All  of  the 
original  test  log  sheets  are  contained  in  the  Appendix. 

The  test  rig  consists  of  a  short  hollow  shaft  supported  by  two  75  mm 
bore  slave  roller  bearings.  The  dual  diameter  roller  bearing,  140  mm  bore, 
is  supported  midway  between  the  two  slave  bearings.  Radial  load  is  applied 
by  pulling  upward  on  the  test  bearing  housing  through  a  large  ring  as  shown 
on  Figure  30.  Oil  to  the  dual  diameter  bearing  is  introduced  through  three 
orifices,  location  No.  1,  and  jeted  to  a  scoop  on  the  test  shaft.  The  oil 
then  travels  to  the  recess  under  the  dual  diameter  bearing,  up  through  the 
bearing,  out  into  the  rig  cavity,  and  exits  the  drain  through  points  2  and 
4.  The  test  shaft  is  driven  by  a  splined  quill  shaft  between  the  rig  and 
the  speed  increaser  gear  box.  The  area  for  the  quill  shaft  drive  is  shown 
with  the  cover  removed  in  Figure  28.  The  entire  test  head  shown  in  Figure  30 
and  a  transition  cone  covering  the  quill  drive  shaft  are  cantilevered  from 
the  speed  increaser  gear  box,  as  can  be  seen  in  Figures  28  and  29. 

Two  prototype  dual  diameter  roller  bearings  were  tested  under  the 
current  contract.  The  first  bearing  was  operated  with  MIL-L-7808  oil 
at  approximately  250°F  inlet  temperature  to  gain  a  feel  for  the  power  losses 
and  general  operation  of  the  bearing  with  a  known  lubricant.  The  first  test 
bearing  was  then  operated  with  Polyphenyl  Ether  5P4E  lubricant  at  approximately 
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500°F  oil  inlet  temperatures  up  to  and  including  3.25  x  10^  DN  operation. 

The  second  test  bearing  was  operated  with  Polyphenyl  Ether  5P4E  lubricant 
at  approximately  500°F  oil  inlet  temperature  and  successfully  completed  30 

g 

minutes  of  continuous  operation  at  3.5  x  10  DN  (25,000  rpm)  with  stabilized 
temperatures.  The  second  test  bearing  »  as  then  subjected  to  a  series  c  ' 
rapid  acceleration/deceleration  tests  and  removed  for  visual  inspection. 

2.  M1L-L-7808,  TEST  BEARING  A1 

A  data  s^r.r_ary  sheet  for  the  operation  of  the  first  test  bearing  s/n 
No.  A1  is  shown  in  Table  XIX.  The  bearing  outer  race  and  front  and  rear 
drain  oil  temperature  increases  as  a  function  of  bearing  speed  for  the 
testing  shown  in  Figure  31. 

Twenty-five  minutes  of  successful  operation  at  3.50  x  10^  DN  was 
obtained  (see  runs  49  and  50),  The  bearing  was  removed  from  the  test  cell 
and  visually  inspected.  The  inner  race,  rollers  and  outer  race  showed  some 
discoloration  due  to  operating  at  300tfF  average  outer  race  temperatures. 

The  rolling  surfaces  looked  to  be  in  ex*  llent  condition  with  only  the 
normal  amount  of  burnishing  to  be  expected  after  completion  of  13  hours 
and  40  minutes  accumulated  running  time. 

The  cage  was  found  to  be  badly  worn  on  the  inner  land  riding  diameters 
over  approximately  a  35  degree  arc  on  both  lands.  Tne  worn  areas  had 
penetrated  the  silver  plating  end  significantly  oatormed  the  steel.  The 
worn  areas  were  directly  under  one  of  tne  cage  pins  used  for  cage  speed 
instrumentation.  It  appeared  that  the  cage  had  excessive  mass  unbalance 
in  operation  which  led  to  the  localized  wear.<g  r  the  inner  land  surfaces. 
Tne  s/n  Nw.  A1  cage  was  not  reused.  The  inner  race,  outer  race  and  roller 
complement  however  were  in  excellent  condition  and  were  reused  .or  the 
second  set  of  tests. 

3.  5P4E,  TEST  BEARING  A1 

s/n  No.  A1  bearing  with  the  s/n  No.  B2  cage  was  assembled  in  the  tost 
rig  for  the  high  temperature  Polyphenyl  Ether  test.  A  data  summary  sheet 
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DATA  SUMMARY  SHEET  FIRST  TEST  BEARING  -  MIL-t-7808  OIL 
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DUAL  DIAMETER  ROLLER  BEARING 
PART  NO.  C3132  S/N  AOOl 


TEST  CONDITIONS: 

LUBRICANT  -  MIL-L-7808G 
LUBE  FLOW  RATE  -  8.5  LB/M  IN 
LUBE  SUPPLY  TEMP  -  250eF 
RADIAL  LOAD  -  480  LBS 
LEGEND: 

Q BEARING  OUTER  RACE 
□  LUBE  OUT  FRONT 
A  LUBE  OUT  REAR 
—  LEAST  SQUARES  LINE 


J _ 1 _  1  . J _ I _ I _ L 

0.5  1.0  1.5  2.0  2.5  5.0  3.! 

BLARING  SPEED -MILLION  i)N 

Figure  21.  Test  Bearing  No.  1,  aT 
vs.  Speed,  MIL-L-7808  G 


is  given  in  Table  XX.  The  temperature  differentials  in  the  bearing  for 
the  outer  race  and  drain  temperatures  are  shown  on  Figure  32.  The  bearing 
was  initially  operated  up  to  and  including  3.0  x  10^  DN  operation.  At 
this  point  the  Polyphenyl  Ether  lubricant  was  being  pumped  back  through  the 
quill  shaft  Figure  30,  to  the  gear  box  and  lixing  with  MIL-L-7808  oil  in 
the  gear  box  system.  The  rig  was  shut  down  and  lube  system  changed  so  that 
the  gear  box  as  well  as  the  test  head  would  operate  on  5P4E  lubricant  and 
would  avoid  the  loss  of  excessive  amounts  of  the  test  oil.  Testing  was 
then  resumed  and  resulted  in  the  second  set  of  curves  as  shown  on  Figure  32. 

Successful  operation  of  the  test  bearing  was  accomplished  up  to  and 
including  3.23  x  10^  DN  operation.  After  approximately  10  minutes  operation 
at  3.50  x  10^  DN  operation  (run  No.  86),  the  bearing  failed  abruptly  as 
evidenced  by  a  seizure  and  a  dead  stop  of  the  test  shaft.  The  failed  test 
bearing  was  then  removed  for  visual  examination.  The  inner  race  had 
separated  and  formed  into  two  circumferential  halves  approximately  in  the 
center  of  the  race.  The  rings  had  also  snapped  open  at  a  location  corre¬ 
sponding  to  one  of  the  key  slots  as  shown  in  Figure  33.  A  considerable 
amount  of  debris  was  trapped  between  the  inner  race  cavity  and  the  shaft. 
Figure  34  is  a  photograph  of  the  rollers,  cage  and  outer  race  assembly 
still  in  the  outer  race  loading  ring.  The  rollers  all  had  a  flat  surface 
on  the  major  diameter  which  presumably  was  caused  when  the  rollers  locked 
up  and  braked  the  inner  race  to  a  complete  halt.  Several  of  the  rollers 
are  show.'  on  Figure  35.  All  of  the  rollers  had  at  least  one  of  the  small 
diameters  separated  from  the  body  of  the  roller  as  shown  in  Figure  35  and 
in  several  cases  both  minor  diameters  had  been  separated.  Subsequent 
metallurgical  examination  of  the  fractured  surfaces  supports  the  fact  that 
the  small  diameters  were  twisted  or  torqued  cut  of  the  major  diameter. 

They  were  not  reversed  bending  fatigue  failures.  Complete  seizure  of  the 
inner  race  creating  sufficient  friction  to  grind  flat  from  the  rollers 
(as  shown  in  Figure  35)  would  create  tremendous  torques  on  the  minor  roller 
diameters.  The  outer  race  is  shown  in  Figure  36  and  the  cage  is  shown  in 
Figure  3?,  The  cage  suffered  serious  distortion  as  would  be  anticipated. 
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SUMMARY  SHEET  FIRST  TEST  BEARING  -  POLYPHENYL  ETHER  5P4E 


DUAL  DIAMETER  ROLLER  BEARINGS 

PART  NO.  C3132  S/N  AOO!  with  B002  CAGE 

TEST  CONDITION: 

LUBRICANT  -  POLYPHENYL  ETHER  5P^E 
LUBE  FLOW  RATE  -  9.5  LB/MI N 
LUBE  SUPPLY  TEMP  -  4808F 
RADIAL  LOAD  -  470  LBS 


LEGEND: 

O  BEARING  OUTER  RACE 
□  LUBE  OUT  FRONT 
A  LUBE  OUT  REAR 
—  LEAST  SQUARES  LINE 
OPEN  1st  RUN  11/6/72 
SHADED  2nd  RUN  11/7/72 


□  ^  C 

o  i 


Figure  34. 


Test  Bearing  No.  1  Outer  Race 
Cage  Roller  Assembly  Failure 


I 


The  test  bearing  A1  accumulated  a  total  of  24  hours  operation  at 
various  speeds  of  which  10  hours  and  20  minutes  were  with  the  Polyphenyl 
Ether  lubricant  at  high  temperature.  The  failure  of  the  inner  rcce, 

Figure  33,  can  be  explained  by  two  possible  failure  modes,  the  first 
being  thermal  lockup  of  the  bearing  where  the  growth  of  the  inner  race 
both  thermally  and  due  to  centrifugal  forces  exceeds  the  growth  rate  of 
the  outer  race  resulting  in  internal  diametral  preload.  Roller  bearings 
are  exceptionally  stiff  and  small  amounts  of  diametral  preload  (on  the 
order  of  0.0002  inches  or  more)  can  result  in  very  high  roller  loadings. 
Greatly  increased  roller  loadings  at  the  inner  race  contact  (see  Figure  17) 
would  then  undoubtedly  fail  in  reverse  bending  fatigue.  The  second  possible 
failure  mode  is  that  sufficient  clearance  did  remain  in  the  bearing  and 
that  failure  cf  the  inner  race  was  due  entirely  to  reverse  bending  fatigue. 
The  finite  element  analysis  of  the  inner  ring  as  shown  on  Figure  13  and 
summarized  in  Table  XX  did  not  indicate  that  excessive  stress  variation 
would  be  anticipated  at  3.5  x  10^  DN  operation.  Therefore,  thermal  lockup 
was  suspected.  The  fact  that  the  bearing  did  operate  for  approximately  10 
minutes  at  the  higher  speed  level,  however,  is  not  indicative  of  thermal 
lockup.  Thermal  lockup  occurs  very  rapidly  after  an  Increase  of  speed. 

The  concern  over  thermal  lockup  was  resolved  by  additional  thermal 
analysis  of  the  test  bearing  operating  conditions.  The  initial  thermal 
gradients  in  the  bearing  for  design  purposes  did  not  include  any  effects  of 
external  energy  being  introduced  into  the  system  because  of  dissipation 
in  the  slave  bearing.  Table  XXI  shows  the  results  of  the  thermal  analysis 
with  external  energy  introduced  at  node  21  on  the  shaft  surface  external 
to  "the  test  bearing.  Also  the  effect  of  insulation  around  the  outer  diameter 
or  outer  surface  of  the  test  head  was  introduced  into  this  analysis.  The 
results  of  the  improved  thermal  analysis  for  the  actual  load,  speed  and 
lubricant  supply  conditions  of  the  test  agree  closely  with  the  observed 
temperatures  of  Table  XX.  The  outer  race  temperature  per  analysis  (node 
18),  was  527.86°F  as  compared  to  average  observed  outer  race  temperatures  of 
520°F.  The  total  AT  to  the  oil  by  analysis  in  Table  XXI  was  (node  11) 
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555°  lass  oil  inlet  temperature  of  470°  resulting  in  an  85°  AT  to  the 
oil.  This  compares  very  favorably  to  the  90°  average  of  the  oil  out 
front  and  oil  out  rear  of  Table  XX.  The  calculated  internal  diametral 
clearance  in  the  bearing  under  the  thermal  gradients  shown  in  Table  XXI 
is  0.0051  inches.  This  number  includes  all  thermal  and  centrifugal  growths. 
A  number  of  variations  in  the  assumptions  of  the  thermal  analysis  such  as 
percentage  air  oil  mixture,  etc.  were  tried  and  found  to  have  little  effect 
on  the  thermal  gradients  in  the  bearing.  Therefore  the  conclusion  was 
reached  that  the  bearing  did  not  fail  due  to  thermal  lockup  but  rather 
failed  due  to  reverse  bending  fatigue  of  the  undercut  inner  race.  The 
finite  element  analysis  summarized  in  Table  IX  would  indicate  that  the 
stress  variations  with  roller  loading  should  not  have  exceeded  the  long 
term  life  value  in  reverse  bending  fatigue,  newever,  the  teat  bearing  at 
the  point  of  failure  had  experienced  39.7  x  10^  cycle*,  of  stress  or  roller 
load  variations  between  rollers  and  inner  races  with  both  MIL-L-7808  and 
the  5P4E  lubricant.  The  values  were  essentially  equal  with  each  type  of 
oil.  The  decision  was  made  to  proceed  with  testing  of  the  second  bearing 
on  the  assumption  that  thermal  lockup  was  not  the  failure  mode  and  to 
limit  the  bearing  load  to  nc  more  than  500  pounds  so  as  not  to  significantly 
increase  the  danger  of  failure  due  to  reverse  bending  fatigue. 

4.  5P4E,  TEST  BEARING  C3 

The  second  test  bearing.  Serial  No.  C3,  was  assembled  in  the  rig 
for  high  temperature  testing  with  5P4E  lubricanc.  A  summary  data  sheet 
of  the  test  is  given  in  Table  XXII  and  the  AT's  as  a  function  of  speed 
are  plotted  in  Figure  38.  This  bearing  achieved  30  minutes  of  temperature 
stabilized  operation  at  3.5  x  10^  DN  (25,000  rprc)  under  a  471  lb.  radial 
load.  The  only  variation  in  test  conditions  over  the  first  test  bearing 
was  the  addition  of  an  auxiliary  jet  towards  the  rear  of  the  bearing  at 
the  inner  race  and  cage  land  contact.  This  auxiliary  jet  supplied  approxi¬ 
mately  2  lbs.  per  minute  oil  flew  of  the  total  reported  in  the  test  log. 

The  test  bearing  was  then  subjected  to  four  rapid  acceleration  and 
deceleration  tests  between  1.5  x  10^  DN  and  3.5  x  10^  DN  in  15  seconds  each. 
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TABLE  XXII 

SUMMARY  SHEET,  SECOND  TEST  BEARING  POLYPHENYL  ETHER  5P4E 


SUPPLY 


Figure  38.  Test  Searing  No.  2  AT  vs. 

Speed,  Polypheny1,  Ether,  5P4E 
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The  second  test  bearing  is  shown  in  Figure  39.  All  of  the  race 
and  roller  operating  surfaces  were  in  excellent  condition.  The  cage  was 
in  good  condition  and  showed  even  burnished  wear  areas  in  the  cage  roller 
to  pocket  contacting  surfaces.  With  the  exception  of  the  riding  land, 
the  cage  was  in  excellent  condition.  The.  silver  plate  was  not  worn  through 
at  ar.y  position  and  very  little  discoloration  was  noticed.  The  cage  inner 
land  riding  surfaces  however  were  discolored  over  a  90°  arc  indicating 
insufficient  lubricant  to  the  cage  lands.  The  inner  race  shows  a  uniformly 
burnished  area  in  the  center  of  the  race  where  the  rollers  contacted. 
Distress  is  noted  in  the  cage  land  riding  areas  on  the  inner  race.  The 
outer  race  was  in  excellent  condition.  There  was  some  indication  in  the 
discoloration  due  to  temperature  that  the  rear  of  the  bearing  was  cooler 
than  the  front  face  of  the  bearing.  This  can  be  attributed  to  the  auxiliary 
oil  jet  directed  at  the  rear  face  of  the  bearing. 

The  Polyphenyl  Ether  5P4E  lubricant  used  in  the  test  was  supplied  by 
the  Aero  Propulsion  Laboratory,  Wr i gh t-Pa t te rs on  Air  Force  Base.  Some 
of  the  fluid  was  new  fluid  and  some  of  the  fluid  furnished  was  reconstituted 
fluid.  It  wa3  noted  that  when  new  5P4E  fluid  was  added  to  the  oil  sump, 
that  foaming  of  the  oil  occurred  for  a  period  of  10  to  20  minutes  and  then 
subsided.  This  foaming  did  not  occur  when  reconstituted  5P4E  fluid  was 
added  to  the  sump.  The  foaming  action  did  not  seem  to  have  any  effect  upon 
bearing  operation. 
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SECTION  VII 


CONCLUSIONS 

The  second  dual  diameter  roller  bearing  successfully  operated  for 
30  minutes  of  stabilized  temperature  operation  at  a  shaft  speed  of  25,000 
rpn.  This  is  3.5  x  10^  DN  operation  of  a  140  mm  bore  bearing.  External 
redial  load  was  471  lbs.  and  the  lubricant  was  5P4E  Polyphenyl  Ether  at 
approximately  500°F  oil  inlet  temperature. 

The  failure  of  the  first  test  bearing  was  attributed  to  reverse  bend¬ 
ing  fatigue  of  the  M-50  steel  inner  race  at  high  speed  and  600°F  operating 
temperatures.  Vie  calculated  stress  levels  were  below  known  allowable 
stress  levels  for  icrg  term  operation.  This  suggests  that  additional 
carefully  controlled  reverse  bending  fatigue  data  is  required  for  M-50 
bearing  steels  for  elevated  temperature  operation.  Any  future  tests  of 
the  dual  diameter  should  use  a  solid  inner  race,  eliminating  the  undercut. 

The  lubrication  method  of  introducing  the  oil  to  a  cavity  in  the 
inner  race  and  through  orifices  directly  into  the  bearing  intering  was 
generally  successful,  however  distress  and  wear  was  observed  at  the 
cage  inner  land  riding  surface. 

The  sue  .essful  operation  of  the  second  test  bearing  establishes 
the  validity  of  the  design  approach.  A  complete  systems  analysis  using 
existing  state-of-the-art  elastohydrodynamic  technology,  fluid  mechanics, 
theory  and  thermal  analysis  techniques,  can  be  used  to  design  a  completely 
new  bearing  concept  and  enable  Initial  prototype  testing  to  be  successful 
at  speeds  and  temperatures  beyond  the  realm  of  current  experience. 
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ABSTRACT 


An  interdiciplinary  systems  analysis  is  presented  for  high¬ 
speed  gas  turbine  engine  mainshaft  roller  bearings  which  will  enable 
the  designer  to  meet  the  demands  for  ever  higher  rotative  speeds  and 
operating  temperatures.  The  latest  elastohydrodynamic  experimental 
traction  data  is  included.  Analytical  results  cite  a  need  for  better 
definition  of  the  rolling  fricticn  portion  of  the  total  traction.  A 
fluid  mechanics  model  for  the  detailed  analysis  of  fluid  drags  is  developed 
based  upon  a  turbulent  vortex-dominated  flow  and  includes  the  effect  of 
lubricant  flow  through  the  bearing.  Typical  cage  equilibrium  solutions 
show  evidence  of  possible  unstable  operation.  A  complete  thermal  analysis 
including  dynamic  and  thermal  effects  upo:-  bearing  dimensions  and  resulting 
clearances  is  also  included.  Heat  transfi  t  coefficients  are  given  in 
detail.  Shaft  power  loss  and  cage  slip  predictions  as  a  function  of 
load,  speed,  and  lubricant  supply  correlate  well  with  available  experi¬ 
mental  data. 
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INI  V1UUCTICH 


The  use  of  e las cohydrodynamic  lubrication  theory  for  rolling 

bearing  analysis  was  first  suggested  by  Dowson  and  Higginson^.  The 

first  use  of  digital  computers  for  the  analysis  of  high  speed  roller 

[21 

bearing  performance  was  by  Harris1  .  The  Harris  paper  presented  an 
analytical  method  to  predict  skidding  in  high-speed  roller  bearings.  The 
influence  of  the  amount  of  lubricant  supply  on  cage  and  roller  motion  in 
high  speed  roller  bearings  was  experimentally  described  and  analytically 
discussed  by  Boness^.  A  refinement  was  introduced  by  Poplawski^  with 
the  solution  of  all  of  the  loaded  rollers  under  the  varying  load  pattern 
resulting  from  external  radial  load.  He  also  introduced  a  fluid  "churning" 
loss  and  cage  pilot  surface  friction  formulations. 

All  of  the  previous  analytical  work  used  the  Dowson-Higginson^ 

isothermal  elastohydrodynamic  (EHD) ,  formulation.  A  recent  comprehensive 

review  of  EHD  theory  and  experimental  data  by  McGrew,  et  al^  resulted 

in  selection  and  recommendation  of  the  latest  and  most  applicable  EHD 

[5] 

formulations.  ‘)his  preliminary  EHD  design  manual  also  presents 
computer  subroutines  for  the  calculation  of  oil  film  thickness  including 
thermal  effects  developed  by  Cheng^  and  EHD  traction  subroutines  fitted 
to  available  disk  machine  experimental  data. 

The  present  paper  presents  a  system  analysis  of  a  high-speed 
roller  bearing  under  radial  loading  as  commonly  used  nn  gas  turbine  engine 
mainshafts.  Oil  film  thickness,  traction,  and  hydrodynamic  pressure  forces 
acting  in  an  elastohydrodynamic  roller-race  contact  are  computed  according 
to  the  latest  available  techniques A  complete  fluid  mechanics  model 
is  developed  in  order  to  more  accurately  predict  the  fluid  drag  losses 
acting  upon  roller,  race,  and  cage  moving  surfaces.  Turbulent  vortex- 
dominated  flow  theory  is  used  to  calculate  boundary  layer  stresses.  Finally* 
a  complete  thermal  analysis  is  described. 

-1- 
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Combining  Che  resulting  thermal  deformations  with  deformations  of  the 
bearing  components  due  to  centrifugal  growth  completely  defines  changes 
in  critical  bearing  operating  parameters  at  each  step  of  the  iterative 
solution  process.  The  thickness  of  an  EHD  oil  film  is  determined  by  the 
temperature  and  viscosity  at  the  inlet  to  the  contact.  This  inlet  oil 
temperature  is  very  sensitive  to  the  temperature  gradients  which  exist 
within  the  bearing.  Oil  viscosity  changes  and  dimensional  changes  through¬ 
out  the  bearing  as  a  function  of  temperatures  are  included.  The  thermal 
solution  as  expected  is  dependent  upon  lubricant  flow  rate,  method  of 
application  (under  inner  race,  thorough  bearing,  etc.)  aid  external  heat 
paths, 

LOAD  AND  SPEED  EQUILIBRIUM 

Elastohydrodynamic  Pressure  Forces  and  Tractions 

The  surface  velocities  and  relative  rotational  speeds  of  rollers 

12] 

and  raceways  are  shown  in  Fig.  1  and  follow  the  nomenclature  of  Harris1  . 
Sliding  velocities,  and  are  defined  as  positive  when  the  surface 

velocity  of  the  ra;e  is  greater  than  the  surface  velocity  of  the  roller 
and  is  equal  to  their  difference.  The  rolling  or  entrainment  velocities 
and  are  always  positive  and  defined  as  the  average  of  raceway  and 
roller  s^.race  velocities.  The  sliding  velocities  and  entrainment  velocities 
are  made  dimensionless  by  introduction  of  the  fluid  entrance  viscosity, 

n  ,  reduced  modulus  of  elasticity,  E*,  and  equivalent  cylinder  radius,  R, 

S  [13 

for  each  contact  as  suggested  by  Dowson  . 

Forces  and  moments  acting  on  a  loaded  roller  are  shown  in  Fig.  2 
and  those  acting  on  an  unloaded  roller  are  shown  in  Fig.  3.  The  effect  of 

cage  orbital  speed,  u  ,  is  to  create  a  dynamic  body  force  on  the  roller 

C  f  71 

acting  radially  outwards.  Figs.  2  and  3  use  D’Alembert's  principle1  J  to 

represent,  the  summation  of  this  roller  body  forces  as  an  inertia  force 

CF*  acting  through  the  center  mass  of  the  rollers. 

CF*  Ew  2  (i) 


-2 


Thus  Figs.  2  and  3  are  free  body  diagrams  of  the  roller  in 
static  and  dynamic  steady  state  equilibrium  with  all  static  and  dynamic 
forces  shown. 

Fluid  pressure  forces  Q. .  and  Q  .  tending  to  translate  the 

ij  oj  m 

roller  relative  to  the  raceways  are  described  by  Dcwson,  etal  and 

f  2 1 

approximated  by  Harris1  for  synthetic  lubricants  as: 

« 

Q  -  18.4  U-Y)  °-°-JU  °-7 

I 

Qoj  -  18.4  (1+y)  G"0,3  Uoj0*7  - 

Traction  forces  T. 4  and  T  .  consist  of  fluid  frictional  drag 

forces  as  described  by  elastohydrodynamic  lubrication  theory.  They  con- 

[21 

sist  of  rolling  and  slicing  friction  forces.  Using  the  Harris  approxi- 


mation  for  and  they  are  expressed  as: 

T  ,  *  ~9.2G~03U  0.7  l  E’R  +-^ 
ij  ij  ij  i 

iul 

V13 

P1J 

T  *  --9.2G-0*3  0  .°*7i  ,  E'R  + 
°j  oj  oj  0 

J^oi 

v 

1  «*■>.*«. 

°j 

The  coefficient  of  traction,  f,  for  the  sliding  portion  has 
been  obtained  by  experiment  from  disk. machines.  This  coefficient  has 
been  curve  fitted  to  existing  data  by  McGrew,  etal^  as  a  function  of 
three  dimensionless  parameters  G^,  G2,  G^  and  is  represented  symbolically 
as: 

f (Te)  -  f(g.,  G  ,  G  ,  30°C  -  0.001  (Te-86°F) 

where  Te  is  the  temperature,  °F,  at  the  entrance  to  the  EHD  contact.  The 
three  dimensionless  parameters  (given  in  the  nomenclature)  express  shear 
rate  effects,  thermal  heating  effects,  and  pressure  viscosity  effects, 
respectively.  Graphs  and  computer  subroutines  for  the  traction  coeffi¬ 
cient,  f  ,  are  given  in  reference,  [5]  and  are  not  reoeated  in- this  paper. 


Ic  should  be  noted  that  this  representation  of  experimentally  obtained 
traction  data  is  extrapolated  for  both  high  and  low  sliding  velocities  to  ■ 
make  the  data  useful  over  a  wide  range  for  the  analysis  of  bearing  systems. 
At  smull  values  of  all  curves  of,  f,  have  a  45°  slope.  This  corresponds 
to  the  fact  that  at  small  sliding  speeds  the  frictional  coefficient  varies 
linearly  with  sliding  speed. 

and  P  are  raceway  normal  forces  transmitted  to  the  roller 
by  the  external  bearing  load.  F  ^  is  a  force  resulting  from  contact 
between  the  roller  and  the  guiding  shoulders  of  the  raceway.  is  a 

symbolic  representation  of  all  of  the  fluid  resisting  moments  acting  on  the 
roller.  F^  is  loading  due  to  the  cage.  A  loaded  roller.  Fig.  2,  which 
drives  the  cage  results  in  a  positive  value  of  F^  as  shown.  An  unloaded 
roller  is  driven  by  the  cage  as  shown  in  Fig.  3.  A  retarding  torque  is 
always  present  as  a  result  of  rubbing  coulomb  friction  between  the  roller 
and  cage  pocket. 


Load  Roller  Equilibrium 

Summation  of  forces  in  the  horizontal  and  vertical  directions. 
Fig.  2,  for  the  loaded  roller  must  equal  zero  for  steady  state  translating 
equilibrium. 

Qij  +  Tij  '  Qoj  “  Toj  "  Fcj  ±  Fadd  "  0 


P.,  +  CF  -  yF  .  -  P  .  «  0 
ij  cj  oj 


Summation  of  moments  about  the  roller  axis  of  rotation,  Fig.  2, 
must  equal  zero  for  steady  state  rotary  motion  equilibrium:  (The  term 
Fadd  resu^-*n8  from  roller  end  contact  with  the  guiding  shoulder  is  des¬ 
cribed  later  by  Eq.  (36). 

,  .  .  „  3(jr, 

X  —  +  T  —  -  uFc  —  -  M  4  0  *  I  - i 

Alj  2  +  oj  2  ytcj  2  rj  f  U  rj  3t 

The  equations  of  equilibrium  for  the  unloaded  roller  are  similar 
and  can  easily  be  derived  from  Fig.  3.  (Also  Ref.  (2]). 
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EHD  Oil  Film  Thickness 


The  average  oil  film  thickness  in  the  Hertzian  zone  is  computed 

t-  aa  (-Uft  ™  a  (-  4  nA(*llo  t*w  <1*1  hllQAi*il  f  AV  1  ^nA  A  A  <  <  1 A  <  An  A 


according  to  the  Dowson^  etal*  1  isothermal  theory  for  line  contact  and 
modified  with  the  Cheng^  thermal  correction  factor*  <f>T.  ^  The  film  thick¬ 


ness  is  determined  by  the  lubricant  velocity,  n  at  the  inlet. 

e » 

.  1.6  G°*6  s  °'7r 

9  °*13 


Roller  Load  Distribution 

The  summation  of  all  inner  race  loadings,  P  ,  equals  the  applied 
external  radial  load,  F  .  For  simplicity , the  solution  will  be  given  in  a 
single  degree  of  freedom,  along  X,  Then: 


F  -  E  P.,  cos  *  *  0  -  cLD 

X  Ij  J 


where 


COS  * 


Cage  Equilibrium 


Steady  state  rotation  of  the  cage  requires  that  the  summation 
of  all  fluid  film  forces  and  roller  to  cage  forces  acting  on  the  cage  be 
in  equilibrium.  The  sum  of  the  moments  of  all  of  the  fluid  film  fo-;ces 
acting  on  the  cage  surfaces  and  taken  about  the  center  of  bearing  rotation 
will  be  represented  symbolically  by  and  will  be  considered  to  be  a  drag 
or  retarding  moment  when  expressed  in  the  positive  sense.  An  inner  land 
riding  cage  will  develop  a  driving  torque  in  the  land,  which  helps  to  dri’-e 
the  cage.  An  outer  land  riding  cage  develops  a  drag  torque  in  the  lands. 

An  unloaded  roller  by  definition  has  a  negative  cage  force  F£,  Fig*  3. 

A  loaded  roller  has  either  a  psotive  driving  cage  force,  Fc,  or  a  negative 
drag  force,  per  Eq.  (5)  and  Fig.  2. 


E  * 

—  £ 

2  U 

J=1 


tV 


Mc  *  0 


Method  of  Solution 


A  flow  chart  of  the  method  of  solution  is  given  in  Fig.  4.  The 
solution  consists  of  six  nested  single  degree  of  freedom  iterative  loops. 
Two  convergence  schemes  are  used.  Load  distribution  solutions  are  solved 
using  classical  Newton-Raphson  matrix  inversion  techniques1  Roller 
speed  and  cage  speed  solutions  use  a  modified  time  step  solution  because 
of  the  multiple  valued  traction  function* 

Equations  (7)  and  (11)  are  written  as  inequalities  where 
the  right  hand  side  is  the  product  of  the  polar  moment  of  inertia  and  the 
first  time  derivative  of  the  rotational  speed.  The  desired  steady  state 
solution  will  occur  when  the  right  hand  side  vanishes.  The  magnitude 
and  sign  of  the  right  hand  side  js  used  to  compute  the  increment  of  cage 
or  roller  speed.  Roller  inertias  are  small,  and  extremely  small  time 
steps  would  be  required  for  solution.  A  modified  time  step  in  the  form 
of  cage  or  roller  speed  steps  proportional  to  the  right  hand  side  (or 
moment  unbalance)  are  used  for  the  solution.  Initial  step  size  is  deter¬ 
mined  by  experience  to  minimize  computer  execution  time.  Changes  in  speed 
regieras  (DN  values)  and  change  of  lubricant  usually  requ'  es  an  adjustment 
in  step  si2s. 


Determination  of  Loaded  Rollers 

The  radial  displacement,  6<^,  at  any  location  angle,  <J>j,  is: 

6$  =  Xcosi* .  -  +  x’  , 

J  J  *  Y 

where  x’  is  a  constant  deviation  from  a  true  circular  profile. 

Y 

A  plus  value  of  x*  means  a  radially  outward  deviation  of  the  inner  race 

or  a  radial  inwaru  deviation  of  the  outer  race  at  vocation  $ ,  or  the 

algebraic  sura  of  deviations  of ’both  races  from  a  true  circular  profile. 

[91 

An  elliptical  out-of-round  race  or  any  race  shape  in  general  1  can  be 
described  by  specifying  values  of  x’  for  every  location  angle, 
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(12) 


Determination  of  a  loaded  roller  is  then  made  by; 


^  i  ‘total  pij  *  0  “'L0ABED 

»  ‘total  pij  »  0  L0ADED 


Eq.  (17)  states  that  if  the  total  radial  displacement,  5$4,  is 

greater  than  the  radial  displacement  value,  of  an  unloaded  roller; 

then  the  roller  is  loaded.  If  6$.  >  6_„..  then  P ../l..  must  equal  or 

j  TOTAL  ij  ij 

exceed  1.0  lb /inch.  This  test  is  made  at  each  roller  location  as  shown  by 
iteration  loop  4,  Fig.  4,  for  the  radial  load  solution. 

Solution  of  a  Loaded  Roller 

The  solution  of  a  loaded  roller  according  to  Eq;  (5),  (6),  and 
(7)  is  accomplished  by  means  of  two  nested  iterations  //5-Iiii'er  Race  Roll 
Load,  P^  and  //6-Roll  Speed,  ,  as  shown  in  Fig.  4. 

The  inner  loop  //6  to  determine  roll  speed  w  (by  means  of  Eq. 

(7))  is  solved  using  the  modified  time  step  technique  described  above. 

This  solution  requires  a  known  value  of  the  inner  race  load,  P^. 

The  inner  race  loading,  P^»  a  loaded  roller  is  determined 
by  using  Newton-Raphson  iteration  (loop  5  of  Fig. -4)  and  assuming  an  initial 
value  of  P^  an  then  calculating  corrections  to  it.  An  initial  value  can 
be  closely  determined  from  the  external  radial  load  and  using  Stribeck’s 
constant  for  the  maximum  roller  load. 

The  cage  to  roller  contact  force  F  is  solved  directly  from  Eq. 
(5)  for  a  given  value  of  w  . .  Roller  speed  is  determined  by  iteratively 
solving  Eq.  (7)  with  a  modified  time  step  or  incremental  values  of 
The  solution  of  Eq.  (7)  is  very  sensitive  to  small  changes  to  torques  about 
the  roll  axis. 


The  pFcj  term  of  Eq.  (6)  is  neglected  in  the  solution.  The 
coulomb  coefficient  of  friction  u  is  usually  less  than  0.1.  The  traction 
coefficient  in  the  EHD  contact,  f(Te),  is  always  less  than  0.1  and  usually 


less  than  .01.  Tht's  the  term  Is  always  less  than  IS  and  usually 
less  than  0.1%  of  P  . 

The  relationship  of  radial  displacement,  x,  of  the  inner  ring 
upon  a  loaded  roller  is: 

5oj  +  Slj  -  *oj  -  hij  +  ScFj  ■*  Stj  -  ‘rotj  -  ltcos*j  +  o  -  Cj  (18) 


Where<$„_  and6  refer  to  a  hollow  roller  and  are  omitted  from 

crj  rj 

Eq.  (18)  when  solving  solid  rollers.  <$  is  the  ring  compression  of  a 

[11]  ^ 

hollow  roller  per  Timoshenko  as: 


<5 
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(19) 


(20) 


Subsequent  values  of  are  then: 


6e 


j  / 


6P 


(21) 


Solution  of  Eq.  (18)  using  corrected  values  of  from  Eq.  (21) 
is  rapid  and  uo  convergence  difficulties  have  been  encountered. 

The  load  carried  by  each  loaded  roller  is  determined  by  the 
Radial  load  solution  loop  No.  4  of  Fig.  4.  Jq.  (9)  is  solved  by  Newton- 
Raphson  iteration  using  the  radial  displacement,  x  of  Eq.  (18)  as  the 
variable.  An  initial  value  of  x  is  computed  using  an  assumption  for  the 
value  of  P^j,  Eqs.  (6,  15,  and  16),  assuming  pure  rolling  (  no  slip)  for 
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FLUID  MECHANICS  MODEL 


Flow  models  for  the  (usually)  turbulent  flow  field  are  briefly 
discussed,  and  then  applied  to  calculation  of  fluid  drags  acting  on  roller 
and  cage  surfaces.  The  influence  of  tne  mixed  component  lubricant-air 
fluid  field  is  treated  with  respect  to  definition  of  wall  shear  stresses 
explicit  in  the  calculation  of  fluid  drags.  Lastly,  the  consequences  of 
the  fluid  drags  acting  on  roller  and  cage  surfaces  are  described  in  relation 
to  their  influence  on  bearing  behavior. 

[Al 

At  le~st  one  previous  investigator,  Poplawski  ,  has  included 
the  effects  of  fluid  drags  in  an  analytical  model  of  roller  bearing  per¬ 
formance.  This  model  lumped  the  effects  of  fluid  drags  into  a  "churning 
loss"  and  suggested  the  use  of  an  effective  air-oil  density  to  improve 
the  analytical  description.  The  Poplawski ^model  demonstrated  signifi¬ 
cantly  better  correlation  with  existing  experimental  data  than  did  the 

f  2i 

original  Harris  model.  The  present  work  examines  in  more  detail 
the  influence  of  fluid  drag  torques  acting  on  individual  surfaces,  and 
accounts  for  the  driving  influence  of  such  torques  on  certain  cage  and 
roller  surfaces  as  well  as  the  retarding  influences  of  such  torques. 

Both  the  energy  dissipation  and  the  convective  heat  transfer  within  the 
bearings  are  closely  associated  with  the  viscous  fluid  behavior  of  the 
lubricant.  The  techniques  employed  in  this  analysis  use  a  combination 
of  empiricial  and  theoretical  treatments  to  define  these  effects. 

From  the  standpoint  of  the  dynamic  solution  of  the  interacting 
rollers,  races,  and  cage,  the  most  important  factor  is  calculation  of  the 
friction  drags  on  the  element  surfaces.  In  general,  this  effect  is 
described  by  an  equation  of  the  form: 

T  -  x  Ar  (23) 

w 
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where  T  is 

A  is 

r  is 

t  is 
w 


the  drag  torque  acting  over  the  element  surface, 
the  surface  area, 

the  reference  radius  from  the  center  of  rotation, 
the  wall  shear  stress 


The  description  of  the  wall  shear  stress  is  a  function  of  the 
fluid  properties,  the  motion  of  the  surface  with  respect  to  the  fluid 
body,  and  the  proximity  of  other  surfaces  influenced  by  the  same  fluid. 
Regarding  the  last  mentioned,  close  clearances  between  rotating  members 
of  the  bearing  have  especially  strong  influences. 

Fritz^2jhas  investigated  the  wall  shear  stress  phenomenon  in 

journal  bearings  with  vortex  turbulent  flow.  Examination  of  typical 

roller-cage  clearance  to  roller  radius  ratios,  and  Reynolds  and  Taylor 

numbers  representative  of  1  x  10°  DN  and  higher  bearing  operation  clearly 

place  such  fluid  flows  in  the  turbulent-vortex  regime  as  defined  by 
fl2] 

Fritz  and  illustrated  for  a  typical  case  in  his  Fig.  3.  For  surface 
finishes  and  bearing  dimensions  typical  of  1  x  10^  DN  and  higher  roller 
bearings,  the  flows  representative  of  cage  surfaces  appears  to  be  in  the 
couette  turbulence  regime. 

Wall  shear  stress  for  a  surface  rotating  in  a  viscous  fluid  is 
defined  by  the  relation  : 

tw  =  f (1/2  o  U2)  (24) 


where  U  ®  mass  average  velocity  of  fluid 
P  -  fluid  density 
f  =  friction  factor 

For  this  correlation,  U  =  — ^ 

2 


The  applicable  friction  factors  from  Fritz 
fitted  for  computer  analysis’  are: 


1.3 


0.539474 


and  curve 


<  25) 
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for  the  vortex- turbulent  correxation  and 


t!  -  3.0 
L 


N. 


Re 


2500 


0.85596 


for  the  couette  turbulent  correlation. 
Here, 

ruC 


N 


N, 


Re  v 
rcoC 


Ta 


wC  /jT 

V  V  r 


and 


fL  *  N 


16 


is  the  Reynolds  number 
is  the  Taylor  number 


is  the  laminar  friction  factor 


Re 


Should  calculation  yield  a  Reynolds  number  less  than  2500,  or  a  Taylor 
number  less  than  41,  a  circumstance  virtually  impossible  in  high  speed 
roller  bearings ,  the  appropriate  friction  factor  then  can  be  taken  as 
the  laminar  friction  factor  f^. 

In  addition  to  the  fluid  drag  acting  on  the  cylindrical  sur¬ 
faces,  the  roller  ends  or  cage  sides  provide  an  additional  source  of 
drag.  Per  Refs.*  [13]  and  [14]  for  a  rotating  disk  wetted  on  both  sides 
the  total  moment  resulting  from  fluid  drag  about  the  disk  centerline  is 
given  by : 


where 


C 

n 


„  1  2  5  r 

“t  "  2  p  u  r  Cn 

3,87/(NRe) for  laminar  flow,  NRg  <  300,000 
0.196/(NRe)1/S  for  turbulent  flow,  NRg  >  300,000 


where  here 
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For  hollow  rollers  and  for  the  cage,  account  must  be  taken  of 
the  areas  pertinent  to  the  calculations  by  modifying  the  radius  factor 
of  the  disk  equation  to  the  form: 


where 


5 

r 


out 


) 


(33) 


r  _  is  the  outer  radius 
out 


r.  is  the  inside  radius 
in 

This  is  at  least  trur  to  a  good  approximation  in  the  case  of  the  roller, 
since  examination  of  Jie  integral  from  which  a  closed  form  solution  of  (30) 
is  derived  shows  that  approximately  75%  of  the  total  moment  contribution 

comes  from  the  area  defined  by  the  outer  50%  of  the  roller  radius.  In 
the  case  of  the  cage,  the  "disk"  torque  acting  on  the  cage  side  surfaces 

is  less  justifiable  on  theoretical  grounds,  but  the  calculated  effect  has 
been  shown  to  be  small  with  respect  to  the  cylindrical  surface  contribu¬ 
tion  in  this  instance. 

Basic  to  the  above  calculations  is  the  realization  that  high¬ 
speed  roller  bearings  never  are  completely  flooded  with  lubricant,  and 
seldom  are  more  than  15  to  20%  full  of  oil  within  the  bearing  for  very 
high  speed  operation.  Indeed,  an  excess  of  oil  within  the  bearing  has 

been  conclusively  shown  by  Boness  to  produce  adverse  effects  on  bearing 

-  £3]  , 

operation  .  Motion  of  the  many  moving  parts  within  the  bearing,  in 
particular  the  pumping  action  of  the  rollers,  will  tend  to  induce  an  air 
oil  mixture  within  the  bearing.  Ilueristically,  it  can  be  argued  that  for  pure 
viscosity  effects  in  a  turbulent  regime,  this  phenomenon  will  have  little 
tendency  to  change  the  fluid  behavior  from  that  of  a  flooded  bearing, 
since  in  general  such  effects  are  confined  to  relatively  thin  boundary 
layers  adjacent  to  the  surfaces.  Further,  the  turbulent  mechanisms  of 
momentum  and  energy  transfer  depends  on  physical  displacement  of  fluid 
particles.  However,  the  density  or  inertia  effects  of  the  air  oil 
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mixture  can  be  drastically  influenced  by  comparison  to  that  of  a  flooded 
bearing.  Hence,  it  becomes  necessary  to  modify  the  density  of  the  fluid 
mixture,  choosing  a  volumetric  basis  for  the  calculation: 


p  *  (p  x  oil  volume  +  p  .  x  air  volume) /Total  Volume 
av  oil  air 


x  ,  defined  as  pQil  x  DECFUL  -  DENS 


oil  Total  volume 


(34) 


The  apparent  density  of  the  fluid  is  taken  as  this  in  roller  drag  calcu¬ 
lations.  Recognizing  that  the  land  regions  and  the  cage  surface  on  the 
close  clearance  land  side  will  be  somewhat  restrictive  to  free  flow  of 
oil  out  such  regions,  the  oil  density  in  this  region  is  taken  as: 


p  , ^  x  DECFUL 
PH1GH  *  0.4  +  0.6  DECFUL 


(35) 


The  fluid  density  in  other  regions  is  taken  as? 

p  x  DECFUL2 
oil 

PL0W  “•  0.4  +  0.6  DECFUL  _  Ut>-' 

These  relations  derive  from  flow  continuity  considerations,  assuming  40% 
of  the  total  flow  will  be  on  average  confined  within  the  higher  density 
region.  While  admittedly  arbitrary,  these  assumptions  have  led  to  excellent 
correlation  with  existing  empirical  data. 


Roller  Fluid  Drag  Torque 

Consider  the  drag  torque  acting  on  the  cylindrical  roller  sur¬ 
face,  as  computed  from  equations  (23)  through  (29),  Trcyl*  In  this  case» 
careful  consideration  must  be  given  to  the  calculation  of  the  character¬ 
istic  radial  clearance  required  in  the  calculation.  From  the  close 
clearance  region  between  the  rollers  and  the  cage,  and  the  interrace  cage 
geometry,  it  is  possible  to  define  a  characteristic,  clearance  between 
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the  rollers  and  the  surrounding  elements.  Figure  5  illustrates  the 
geometry.  This  clearance  is  defined  by: 

c  ■  CCH  ■  h  |  c<e>d8  <37) 

Circumference 

Five  distinct  zones  must  be  considered  in  the  calculation  on  each  semi- 
cylindrical  half  of  the  roller.  Rollers  are  assumed  to  be  in  contact  with 
either  the  leading  cage  pocket  surface  (driving  roller)  or  the  lagging 
cage  pocket  surface  (driven  roller).  Generally,  rollers  carrying  load 
are  driving  rollers  and  unloaded  rollers  are  always  cage-driven.  Com¬ 
paring  the  two  roller  surfaces,  therefore,  the  only  zone  not  symmetric 
with  respect  to  the  radial  ray  through  the  roller  center  is  Zone  III, 
where  the  entire  diametral  pock°c  clearance  is  taken  as  10%  of  the  roller 
radius.  Calculation  of  the  clearance  is  then  straightforward,  given 
dimensions  of  the  cage,  rollers,  and  races. 

Similarly,  the  drag  torque  acting  on  the  roller  ends,  T 

rend 

is  computed  directed  from  Equations  (30)  through  (33).  For  all  roller 
calculations,  the  appropriate  density  is  computed  from  Equation  (34). 

While  not  totally  fluid  mechanical  in  origin,  the  sensitive 
nature  of  the  roller  equilibrium  solution  requires  consideration  of  the 
influence  of  the  contact  betv.en  the  roller  and  the  guiding  shoulders 
of  the  cage  in  the  torque  equilibrium  sr  ition.  Computer  results  show 
this  effect  to  contribute  generally  5  to  10%  of  total  roller  drag  moment. 
Further,  for  a  roller  guided  on  the  outer  race  shoulder,  this  moment  can 
either  drive  or  retard  the  roller,  depending  on  the  distribution  of  relative 
local  velocities  between  roller  and  shoulder.  The  calculation  requires 
the  assumption  of  a  coulomb  coefficient  of  friction  between  the  roller 
and  shoulders,  generally  taken  as  0,05  to  0.08  for  a  weJl  lubricated 
contact.  While  strictly  an  integral  over  the  contact  area,  the  moment 
can  be  closely  approximated  by  dividing  the  contact  area  into  a  number 
of  horizontal  and  vertical  strips  or  lamina  as  shown  in  Figure  6  for  a 
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roller  guided  on  the  outer  race.  A  similar  situation  exists  for  a  roller 
guided  on  the  inner  race,  but  the  appropriate  relative  kinematics  must 
be  used.  Further,  for  the  inner  race  guided  roller,  the  local  race 
speed  always  exceeds  the  local  roller  speed,  and  hence  always  produces 


a  net  driving  '  ■  jue  on  the  roller.  The  moment  contribution, 
the  contact  tc  t  roller  moment  is  then  given  by 

f  1 


I  „ .  u  J 

add  c  A_ 
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1  (V  -  V 
N=»l  v  Ri  Hji 


V  -  V 
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S^2 

rvJk 


of 


(38) 


The  second  term  always  produces  a  retarding  moment,  while  the  sign  of 
the  first  term,  at  least  for  an  outer  race  guided  roller,  is  dependent 
on  the  velocity  distribution  over  the  contact  region. 


The  resultant  vertical  force  is  self-equilibrating,  since  the 
velocity  distribution  is  asymmetric  with  respect  to  the  radial  ray  through 
the  roller  center.  However,  the  horizontal  resultant  can  be  either  posi¬ 
tive  or  negative,  and  contributes  to  the  contact  force  between  the  roller 
and  the  cage,  as  evidenced  in  Equation  (5).  This  resultant  force  is  given 
by 


add 


PcFJGf 


|VRi  "  W 
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Vjd 


(39) 


Note  that  these  quantities  are  functions  of  velocity,  and  must  be  deter¬ 
mined  at  each  roller  iteration.  Since  is  constant  in  any  given  pro¬ 
blem,  the  derivatives  of  and  T_^  are  not  involved  in  any  of  the 

Newton-Raphson  i terations . 


The  contact  force  between  the  roller  ends  and  the  guiding  race 
shoulder,  Fjc  acts  in  the  direction  of  the  roller  axis.  Such  a  force  can 
arise  from  a  variety  of  sources.  Among  these  are  differential  pressure 
acting  across  the  bearing,  cage  angular  misalignments,  roller  skewing 
due  to  gyroscopic  moments,  and  unbalances  in  the  various  rotating  mem¬ 
bers  which  give  rise  to  dynamic  forces  between  the  roller  and  the  guiding 
shoulder. 
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In  this  formulation  of  the  problem,  the  contact  force  was  assumed 
to  arise  from  c'  rential  pressure  resulting  from  crossflow  of  lubricant 
through  and  oui  the  bearing.  This  axial  crossflow  was  assumed  to  be 
superimposed  on  L..e  swirling  turbulent  flow  pattern  existing  in  the  annular 
spaces  between  the  rollers,  cage,  and  races.  In  traversing  the  distance 
auwss  the  bearing,  the  resistance  to  flow  in  the  roller  axis  direction, 
and  hence  the  differential  pressure  across  the  bearing,  will  depend  not 
only  on  the  axial  flow,  but  also  on  the  flow  in  the  tangential  direction. 

We  assume  the  same  turbulent  flow  mechanism  to  be  operating  in  both  the 
axial  and  tangential  directions.  This  leads,  after  some  manipulation, 
to  a  differential  pressure  relation 


fp  £  Eu  LV 
P  — c 
16  A 


(40) 


where 


£p  is  the  peripheral  length  of  a  typical  crossflow  passage 


E  is  the  bearing  pitch  diameter 

V  is  a  "length  of  peripheral  surface  averaged"  circumferential 
velocity 

is  the  cage  angular  speed 
L  is  the  bearing  length 

Ax  is  the  area  of  a  typical  crossflow  passage 

[13] 

Calculation  of  the  friction  factor  uses  the  Blasius  relation 


where  here 


f  =  0.3164/(^/n^') 


(41) 


It  is  assumed  that  this  differential  pressure  acts  uniformly  on  the  cage 
side  and  the  roller  disk  area  exposed  to  the  crossflow.  Defining  this 
area  as  Ap,  finally, 


T3a  ■  !>>  \ 


which  is  the  force  appropriate  to  Equations  (38)  and  (39). 


(43) 


The  total  fluid  drag  moment  ,Mr  ^  acting  on  the  roller  in  Equation 
(7)  is  then  calculated  as: 

M  ,  =  T  ,  +  T  .  +  T  . .  (44) 

rj  reyl  rend  —  aad 


Note,  while  the  first  two  terms  are  always  positive  in  this  formulation, 
the  algebraic  sigh  of  T  ^  ii  determined  by  the  relative  velocity  dis¬ 
tribution,  as  discussed  above. 


Ci-je  fluid  Drag  Torque 

The  fluid  drag  torque  acting  on  the  cage  is  composed  of  the 
fluid  torques  acting  on  the  inner  and  outer  cage  surfaces,  the  cage  sides, 
and  the  cage  land  surface.  Its  influence  in  the  equilibrium  solution 
for  the  cage  is  indicated  in  Equation  (11)  by  the  quantity  Hc« 


Calculation  of  the  four  component  drag  torques  proceeds  directly 
from  Equations  (23)  through  (29)  for  the  inner  and  outer  cage  surfaces 
and  for  lands,  and  from  Equations  (30)  through  (33)  for  the  cage  sides. 

The  appropriate  radii  and  clearances  for  these  calculations  are  the  local 
cage  surface  radii  and  clearances  between  adjacent  surfaces.  Fluid 
densities  in  the  various  regions  are  computed  from  Equation  (35)  and 
(36).  Dimensional  changes  with  speed  and  thermal  effects  resulting  from 
dissipation  are  insignificant  in  these  calculations,  except  in  the  case 
of  close  land  clearances. 

The  total  drag  moment  acting  on  the  cage  is  then  computed  from 

M  *  T  -  T  .  +T  ..  +T,  .  (45) 

c  cout  cin  cside  —  eland 
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Note  that  since  the  total  drag  moment  positive  sense  is  in  the  direction 
of  retarding  cage  motion,  the  outer  surface  and  side  torque  terms  are 
always  positive.  The  inner  surface  drag  torque  always  acts  to  increase 
cage  speed  and  is  hence  always  negative.  The  sign  of  the  land  surface 
torque  is  positive  for  an  outer  land  riding  cage  and  negative  for  an 
inner  land  riding  cage. 

Influence  of  Fluid  Dra?  Torques  on  Solution 

Figure  7  compares  the  predictions  of  three  analytical  models 

with  samples  of  previously  published  test  data  of  rolle."  bearing  cage 

slip  as  a  function  of  applied  radial  load.  The  present  model  demonstrates 

almost  exact  correlation  with  this  data  across  the  load  range  for  which 

data  is  available.  Close  correlation  is  shown  by  the  Poplawski ^ model, 

f  91 

which  uses  a  simplified  fluid  drag  correlation.  The  Harris  model 
matches  closely  the  trend  of  the  test  data,  but  differs  significantly  in 
the  magnitude  of  slip  predicted  across  the  load  range.  In  addition  to 

the  correlation  shown  for  an  outer  land  riding  cage,  the  present  model 

* 

has  also  accurately  predicted  the  experimentally  measured  cage  slip  of 
an  inner  land  riding' cage  bearing. 

Power  dissipations  predicted  by  the  model  closely  match  available 

** 

rig  data.  Shaft  horsepower  dissipation  is  computed  by  summing  all 
torques  acting  on  the  inner  race.  The  prediction  and  distribution  of 
these  dissipative  losses  is  essential  for  accurate  determination  of 
thermal  gradients  and  operating  temperatures . -  Predictions  for  two  con¬ 
ventional  roller  bearings  analyzed  for  comparison  are  summarized  below. 
Operation  was  approximately  2.0  x  10^  DM. 


Experimental  Power 
Dissipation,  Hp 


Computer  Power 
Dissipation,  Hp 


130  MM  Bore  5.80  6.94 

115  MM  Bore  4.25  3.89 


* 

Data  courtesy  of  Boeing  Co.,  Vertol  Division,  Philadelphia,  Pa. 

■hit 

Data  courtesy  of  Midwest  Aero  Industries,  Div.  Pure  Carbon  Co.,  St.  ha<ys, 
Pa.  -20- 
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Both  bearings  were  operated  at  near  synchronous  condition. 

This  calculation  is  extremely  sensitive  to  the  amount  of  lubricant  within 
the  bearing,  estimated  at  between  10  to  15%  by  volume  by  the  test  agency. 

The  100  MM  bearing  was  analyzed  for  15%  oil-to-air  ratio  within  the  bearing, 
and  the  115  MM  bearing  was  analyzed  for  a  10%  oil-to-air  ratio.  It  will 
be  noted  that  the  results  bracket  the  experimental  data,  being  high  for 
the  higher  oil/air  ratio  and  low  for  the  lower  oil/air  ratio.  A  rppeat 
of  this  analysis  taking  both  bearings  at  about  a  13%  oil/air  ratio  would 
yield  results  matching  the  experimental  data. 

The  foregoing  illustrates  generally  the  sensitivity  of  the 
solution  to  fluid  drag  influences.  As  further  illustration,  consider 
FigureB  which  compares  two  cage  error  curves  as  functions  of  cage  speed 
for  all  other  conditions  the.  same.  For  this  relatively  low  speed  applica¬ 
tion,  a  difference  of  10%  in  oil/air  ratio  produced  a  2%  difference  in 
cage  slip, 

The  shape  of  the  cage  error  curve,  of  Figure  8,  provides  a 

possible  further  insight  into  the  bi-stable  operation  of  bearings  noted 

f  151 

by  some  investigators.  Dr.  Ehrich  J  noted  this  phenomenon,  and 
hypothesized  its  origin  as  systematic  changes  in  bearing  clearances 
coupled  with  dynamic  phenomenon  of  shaft  rotation.  The  cage  error  vs 
speed  curve  always  has  the  same  characteristic  shape,  for  which  three 
solutions  are  possible  when  slip  occurs.  One  root  is  always  close  to 
synchronous,  where  the  steep  slope  of  the  error-speed  curve  requires 
successive  predictions  at  small  slip  valves  to  obtain  a  solutions. 

Other  solutions  are  possible  if  the  error  curve  changes  sign  at  less 
than  synchronous  speed  (two  possible  additional  solutions),  or  if 
the  error  curve  is  tangent  to  the  abscissa  of  one  point.  Since  the  magnitudes 
of  the  errors  are  strong  functions  of  speed  and  load,  and  to  a  lesser  extent 
of  oil/air  ratio,  it  is  possible  that  for  some  combination  of  these 
variables  the  bearing  operation  may  oscillate  between  two  closely  spaced 
roots.  This  might  occur  as  thermally  induced  dimensional  changes  in 
internal  diametral  and  close  land  clearances  change  the  operating 
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conditions  of  the  bearing.  Juch  operation  has  been  demonstrated  using 
the  present  model  with  thermal  and  dynamic  effects  considered.  One 
computer  run  found  about  3%  slip  on  the  first  iteration  (without  thermal 
gradients) ,  but  sought  the  synchronous  solution  when  the  results  of 
thermal  gradients  resulting  from  the  first  iteration  were  reflected  in 
the  second.  While  not  conclusive  at  this  time,  this  subject  bears 
extensive  further  investigation,  with  comparison  to  actual  test  data  as 
it  becomes  available. 


DIMENSIONAL  MODEL 


The  model  accounts  for  changes  which  occur  during  bearing 
operation  as  a  result  of  dynamic  (centrifugal  effects)  and  thermal 
deformations.  The  effects  considered  to  be  most  significant  are: 


•  Changes  in  bearing  internal  diametral  clearance  (IDC)  with 
dynamic  and  thermal  deformations  and  with  shrink  fit  between 
the  inner  race  and  the  shaft. 

•  Changes  in  bearing  internal  diametral  clearance  (IDC)  due  to 
changes  in  roller  (especially  hollow  rollers)  diametor  with 
dynamic  and  thermal  growth. 

•  Changes  in  close  clearance  guiding  cage/race  lands  with 
thermal  and  dynamic  deformations,  and  for  an  inner  land 
riding  cage,  changes  in  the  clearance  with  shrink  fit 
between  shaft  and  inner  race. 

The  uae  of  simple  closed  form  expressions  to  predict  the 
dynamic  and  thermal  dimensional  changes  occurring  in  bearing  members  is 
desirable,  since  this  greatly  simplifies  the  system  analysis.  Predictions 
of  such  deformations,  on  the  other  hand,  must  be  more  than  order  of  mag¬ 
nitude  accurate  to  realistically  describe  the  influence  of  such  operation 
dependent  variables  on  bearing  performance.  In  the  process  of  model 
selection,  estimates  of  deformations  from  available  closed  form  solutions 
were  compared  to  the  results  more  sophisticated  analyses,  using  finite 
element  computer  solutions  of  typical  bearing  element  geometries. 

Dynamic  Deformations 

The  dynamic  deformation  models  selected  for  the  races,  rollers, 

and  cage  are  given  in  Table  l  and  are  obtained  from  standard  theory  of 

elasticity  sources. ^  an<*  ^  Comparison  of  finite  element  solutions 

for  a  0.50  ’’diameter  hollow  roller  with  a  0.25’’  diameter  hole  showed 

closest  agreement  with  the  generalized  plane  strain  solution.  At  roller 

speeds  of  150,000  RPM  and  higher,  a  ±42  bulging  and  contraction  of  the 

outer  surface  near  the  roller  edges  was  calculated  in  the  finite  element 

solution.  This  cannot  be  predicted  by  any  of  the  closed  form  solutions, 

however,  since  it  is  common  praccice  to  crown  rollers,  this  variation  is 
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not  considered  significant.  Fcr  a  solid  roller,  little  comparative 
difference  was  noted  for  any  of  the  closed  form  deformation  relations, 
and  hence  a  plane  stress  relation  was  used  for  simplicity. 

Both  the  inner  and  outer  surface  deformations  of  the  inner 
race,  are  significant.  The  deformation  of  the  inner  surface  contributes 
to  change  in  the  shrink  fit  between  the  race  and  the  snaft.,  and  deformation 
of  the  outer  surface  influences  both  the  bearing  IDC  and  the  land  clearance 
of  an  inner  land  riding  cage.  The  closed  form  dynamic  deformation  relation 
best  describing  the  behavior  of  both  race  surfaces  was  found  to  be  the 
plane  strain  relation.  These  relations  generally  predict  deformation  within 
5%  of  the  finite  element  results  even  at  extreme  speed,  providing  race 
undercuts  or  other  irregularities  such  as  guiding  shoulders  are  small. 

The  shaft  outer  surface  deformation  is  also  significant  in  that 
it  contributes  to  change  in  shrink  fit.  While  not  checked  against  a 
finite  element  model,  the  shaft  conforms  to  the  basic  assumptions  of  the 
plane  strain  relation. 

The  deformations  of  the  cage  land  surfaces  are  significant  iu 
the  determination  of  clearances  required  in  the  calculation  of  cage  land 
fluid  viscous  drag.  Agreement  of  closed  form  predictions  with  the  results 
of  finite  element  modeling  were  within  15  to  20%  for  this  deformation  for 
all  closed  fora  solutions.  However,  the  calculation  of  the  fluid  viscous 
drag  has  been  determined  to  be  insensitive  to  this  variation,  changing 
only  about  2%  at  maximum  speed.  This  larger  deviation  is  the  result  of 
the  inability  of  the  closed  form  solutions  to  account  for  the  difference 
in  cage  surface  and  land  diameters  or  to  account  for  bending  deformations 
resulting  from  the  additional  mass  of  the  land  protrusions.  Over  the 
relatively  narrow  land  width,  the  land  surface  deformation  is  fortunately 
quite  uniform,  as  predicted  by  finite  element  analysis.  The  closed  form 
solution  best  matching  the  finite  element  results,  within  the  iimitations 
discussed,  was  the  plane  stress  formation. 

Shrink  fit  (A)  between  the  inner  race  and  shaft  act3  to  de¬ 
crease  the  bearing  internal  diametral  clearance  (IDC),  as  compared  to 
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the  IDC  of  an  unmounted  assembled  bearing.  Also,  for  an  inner  land 
riding  cage  bearing,  the  shrink  fit  decreases  the  close  land  clearance. 
An  extensive  treatment  of  .this  subject  is  contained  in  references  [18] . 


Thermal  Deformations 

The  model  used  for  calculations  of  thermal  deformations 
assumes  an  axisymmetric  temperature  distribution,  uniform  end  deforma¬ 
tions  and  cylindrical  geometry.  The  radial  deformation  of  the  ith  bear¬ 
ing  surface,  having  inner  radius,  r^,  and  outer  radius,  rQ,  is  given  by: 
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The  above  relations  are  known  to  be  quite  accurate  within  the 
governing  assumptions-.  No  attempt  has  been  made  to  date  to  confirm  ir. 
detail  cheir  accuracy  for  this  application.  Comparison  of  thermal  growth 
predictions  to  the  limited  empirical  data  available  to  the  authors  shows 
reasonable  agreement  of  the  magnitude  of  these  deformations. 

When  applying  the  above  to  prediction  of  the  thermal  growth  of 
the  outer  race,  it  must  be  noted  that  no  radial  constraint  is  included 
In  the  formulation.  This  is  typical  of  the  test  environment  of  roller 
bearings  where  loading  is  accomplished  through  an  external  yoke,  but 
growth  of  the  cuter  race  is  otherwise  unconstrained.  Other  application 
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of  this  type  of  analysis  must  consider  any  radial  constraint  appropriate 
to  the  application,  particularly  because  of  the  possibility  of  bearing 
thermal  lock-up. 

From  the  above,  it  follows  that  when  the  surface  temperatures 
of  all  bearing  elements  are  known,  the  thermal  deformations  of  the  bearing 
members  are  completely  defined. 

THERMAL  ANALYSIS 
Problem  Formulation 

The  temperature  field  within  the  bearing  must  be  known  in  order 
to  compute  thermal  deformations.  With  the  basic  assumption  of  axisyrametric 
temperature  distribution,  the  techniques  of  performing  such  an  analysis 
are  well  known.  Harris  [19]  has  applied  heat  transfer  techniques  to  the 
determination  of  bearing  temperature  distributions.  The  technique 
presented  here  will  basical  .y  follow  this  presentation,  with  additional 
emphasis  on  the  calculation  and  distribution  of  power  dissipation  within 
the  bearing  and  on  calculation  of  convective  conductances  or  heat  transfer 
coefficients  appropriate  to  the  calculation. 

All  three  inodes  of  heat  transfer,  conduction,  convection,  and 
radiation  can  operate  between  bearing  members  and  between  the  bearing 
and  its  environment.  Conduction  is  a  linear  phenomenon.  Convection  is 
basically  non-liniar  in  that  the  convective  conductance  is  strongly 
velocity  and  temperature  dependent.  However,  in  the  iterative  solution 
technique  used,  Fig,  4,  velocity  is  fixed  at  each  step  of  the  iteration. 
Temperature  dependence  is  mainly  reflected  in  fluid  property  changes. 

Since  the  iterative  scheme  employs  the  updated  temperatures  from  the 
previous  temperature  iteration  soli  .ion  in  the  convective  conductance 
calculations  for  the  current  temperature  solution,  the  convective  mode 
formulation  is  linear  at  each  temperature  iteration.  Radiation  heat 
transfer  depends  on  differences  of  forth  powers  of  surface  temperatures, 
and  is  hence  always  non-linear.  The  forms  of  the  conductive,  convective. 
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and  radiation  heat  transfer  equations  respectively  are: 


f  2nkl 
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AX.  ;  hollow  cylindrical  body  of  length 
^  1  in  radial  direction 
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where 


Boltzmann  radiation  constant 


th 


4  and  £  are  the  emissivities  associated  with  the  a  and 
m  n  th 

n  surfaces,  respectively. 


F 


m,  n 


**  view  factor  of  nC^  surface  seen  from  m^  surface 


These  equations  are  basic,  and  are  derived  in  any  standard  heat  transfer 
text,  such  as  reference  [20] . 

The  radiation  mode  of  heat  transfer  becomes  significant  only 
for  temperature  differ  >  _s  of  the  order  of  100°F  and  higher.  For  normal 
bearing  operation,  such  conditions  do  not  exist  within  the  bearing.  How¬ 
ever,  heat  transfer  by  this  mode  can  be  significant  between  the  bearing 
and  its  environment,  and  must  be  included. 

Formulation  of  the  heat  transfer  problem  requires  discretizing 
the  bearing  elements  into  a  set  of  nodes,  and  writing  a  Kirchkoff  energy 
balance  at  each  node.  The  power  dissipation  occurring  within  the  bearing 


acts  as  Che  driving  potential,  and  must  be  computed  and  distributed  among 
the  nodes.  Fig.  9  defines  the  nodal  pattern  used  in  this  analysis.  In 
this  model,  provisions  are  made  for  coolant  flow  through  the  outer  race 
and  through  the  shaft,  and  for  lubricant  to  be  flowing  both  through  the 
bearing  and  under  the  inner  race.  Any  of  these  flows  except  the  flow 
through  the  bearing  may  be  set  to  zero.  The  coolant  through  the  outer 
race,  and  the  lubricant  flows  under  the  race  and  through  the  bearing  are 
assumed  to  mix  external  to  the  bearing,  and  to  wash  the  inner  surface  of 
the  chamber,  the  shaft,  and  the  exterior  of  the  bearing.  The  environment 
external  to  the  bearing  assumes  the  assembly  to  be  enclosed  in  a  test 
chamber,  and  to  be  loaded  through  an  external  yoke  extending  through  the 
chamber.  The  governing  equations  are  then  of  the  form,  for  nodes  i  ■  1 
to  21* 


21 


l 

J-l 


(50) 


Here  is  the  energy  transfer  between  nodes  i  and  j,  and  the  are  the 
dissipations,  external  energy  inputs,  energy  transfer  to  coolant  or 
lubricant  flows,  or  losses  from  a  typical  node  i.  is  the  total  number 
of  0^  at  node  i.  Table  2  illustrates  the  flow  of  energy  between  all 
nodes,  and  defines  the  quantities  at  each  node.  The  explicit  forms  of 
the  equations  defined  generally  by  (50)  are  derivable  directly  from 
Table  2.  Note  that  all  energy  fluxes  other  than  losses  are  assumed  to 
enter  each  node.  The  temperature  differences  of  equations  (47),  (48) 
and  (49)  may  thus  be  uniformly  taken  at  (T^  -  T^)  at  each  node  i. 

The  energy  balances  of  equations  (50)  are  equated  to  the  error 
functions  e^,  where  approaches  zero  for  the  steady  state  equilibrium 
set  of  temperatures.  At  each  temperature  iteration,  Fig.  4,  in  the  bearing 
solution,  <^is  non-linear  set  of  equations  is  solved  by  Newton-Raphson 
iteration  techn.'nues.  As  an  initial  guess  of  the  temperature  field  to 
start  the  iterative  process,  setting  all  temperatures  equal  to  the  previous 
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value  of  bearing  oil  temperature,  T ^ ,  has  been  found  to  produce  rapid 
convergence.  To  implement  the  Newton-Raphson  solution  scheme,  the  set  of 
equations  (50),  with  the  H. ^  quantities  expressed  in  terms  of  equations 
(47),  (48),  and  (49)  must  be  rewritten  in  the  form 
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Here  the  A„  are  thermal  admittances  for  terms  originating  from  convection 
and  conduction,  and  A*^  are  thermal  admittances  for  terms  originating 
from  radiation.  S|  is  the  total  number  of  quantities,  at  node  i  not 
explicitly  containing  temperature,  and  are  constants  originating  from 
such  sources  as  products  of  lubricant  mass  flow  rates,  specific  heats, 
and  inlet  temperatures.  The  non-linear  set  of  equations  (51)  for  i  *  1 
to  21  is  then  directly  solvable  for  the  temperature  field  using  the  Newton- 
Raphson  solution  technique. 


Power  Dissipation 

As  previously  stated,  the  driving  potentials  which  result  in 
thermal  deformations  are  the  power  dissipations  occurring  within  the 
bearing.  These  arise  from  fluid  boundary  layer  and  coulomb  friction 
phenomenon.  The  distribution  of  the  component  dissipations  is  shown  fn 
Table  2.  The  total  power  expended  within  the  bearing  must  appear  at  the 
shaft  surface,  since  this  is  the  driven  member,  and  is  easily  calculated 
as  the  sum  of  all  roller  to  inner  race  tractions  and  the  sum  of  the  fluid 
drags  acting  on  the  inner  race  and  shaft  surfaces  (including  land  torques 
for  an  inner  land  riding  cage).  Shaft  power  loss  is  an  important  design 
factor  and  one  that  until  new  could  only  be  estimated  empericaJly  based 
upon  test  experience. 

Power  dissipations  at  specific  nodes,  shown  in  Table  2,  are 
also  easily  calculated.  Roller  to  race  contact  dissipation  is  the  product 
of  the  EHD  traction  and  the  sliding  velocity  in  the  contact.  Fluid  drag 
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moments  are  a  ribed  by  the  fluid  mechanics  model. 

Ext  ’,1  energy  Input  such  as  energy  dissipations  by  slave 
bearings  supporting  the  shaft  in  test  apparatus,  other  supporting  bearings 
in  an  engine,  energy  input  from  turbine  gasses,  etc.,  must  be  calculated 
separately  and  introduced  as  heat  sources  or  sinks  at  the  appropriate 
nodes  of  Table  2. 

Heat  Transfer  Coefficients 

Correct  values  of  the  heat  transfer  coefficients,  or  more 
properly  the  convective  conductances,  used  in  the  thermal  analysis  are 
essential  to  obtaining  a  representative  solution.  While  standard  tech¬ 
niques  of  determining  these  quantities  can  be  used  for  non-rotating  sur¬ 
faces,  the  shaft,  cage,  and  roller  surfaces  pose  a  special  problem. 

1211 

3j orklund  has  investigated  the  problem  for  the  case  of  concentric 
rotating  cylinders.  The  techniques  and  correlations  of  this  source  will 
be  used  here. 

For  two  concentric  rotating  cylinders,  with  the  annular  space 
between  them  containing  a  viscous  fluid,  the  critical  Taylor  number  for 
both  is  given  by 


(l-n) 


fl-  n(i+  ~  )T 

-  ' 
0,0571  -^-  0.625  — 

L  r5  J 

l-n  r- 

3  - 

m 

+  .00056 _ 

ir?-  °-625  -r 


(5 


2) 


where 


u 

n  *  Speed  ratio 


(53) 


C  *  radial  clearance  between  inner  race  and  inner  cage  surfaces 
or  between  inner  race  surface  and  cage  land  surtace  for  an 
inner  land  riding  cage. 
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The  Taylor  number  for  this  condition  is  given  by 


NTa  -  /|~ Otwc/B) 


(54) 


where 


6  *  kinematic  viscosity  of  fluid 
R  «  radius  of  inner  rotating  member 

fl  for  inner  race  to  cage  condition, 

to  — « 

wc  for  outer  race  to  cage  condition 
The  Nusselt  number  for  pure  conductance  through  the  fluid  is  given  by 


(55) 


and  the  ratio  of  the  local  Nusselt  number  to  the  conduction  Nusselt  number 
is  given  by 


V  .11 

N  - 

Ta 

flL  -  N_ 

Ta  Ta 

crn  crn«0 

[l  -  3.5  — 1 
r5 

NNu  .  * 

cond 

41.1  + 

[NTa  “  NTa  r  1 

crn  Vo 

,  4 

3-5  h 

‘  4 

1/2 


(56) 


Finally,  the  convective  conductance  is  given  by 


h 


where  k^  *  thermal  conductivity  of  the  fluid. 
For  the  inner  cylinder  rotating  only, 


Nu 


0.175 


Nu 


1/2 


cond 


(57) 


(58) 


applies  as  an  alternate  to  equation  (53). 
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Here  the  appropriate  C  is  the  radial  clearance  between  the  outer  race 
and  cage  outer  surface  or  cage  lands  for  an  outer  land  riding  cage. 
Equations  (52)  through  (58)  can  be  used  directl>  to  obtain  the  appropriate 
coefficients  for  the  race  and  cage  cylindrical  surfaces.  Over  the  ver¬ 
tical  cage  surfaces,  an  average  of  the  coefficients  of  the  bounding 
cylindrical  surfaces  is  recommended.  These  relations  also  define  the 
coefficient  values  on  the  roller  surfaces,  using  equation  (58)  and  the 
average  roller  clearance  defined  by  equation  (37).  In  all  of  the  above 
calculations,  the  fluid  density  used  for  the  calculation  of  the  kinematic 
viscosity  should  be  computed  as  discussed  in  equations  (34-36) .  The 
variation  of  fluid  properties  with  temperature  must  also  be  considered, 
updating  the  fluid  properties  at  each  step  of  the  iterative  solution. 


Coefficients  for  the  outer  bearing  surfaces  and  the  chamber 

inner  and  outer  surfaces  are  computed  from  a  correlation  suggested  by 
[271 

Jakob1  J  for  natural  convection  of  fluids  in  contact  with  horizontal 
cylinders 


where 


0.525 


D;  Yg  At 


0.25 


*»  diameter  of  cylinder 
Y  *  coefficient  of  volumetric  expansion 
g  »  gravitational  acceleration 
*  thermal  diffusivity 

At  ■  absolute  value  of  temperature  difference  between 
surface  and  bulk  fluid  temperature 


(59) 


For  condition  inside  the  chamber,  the  fluid  properties  appropriate  to 
the  calculation  are  those  of  the  lubricant,  since  the  bearing  outer  sur¬ 
face,  the  shaft,  and  the  inner  chamber  wall  are  assumed  to  be  washed 
with  the  mixed  oil.  The  temperature  difference  should  be  taken  between 
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the  surface  temperature  in  question  and  the  mixed  lubricant  temperature, 
which  is  computed  from 


where 


T 

ave 


m. 


pb  X6  *  mll  P11  T11  *  m16  P16  X16 

cp6  +  “h  cPn  +  “i6  Se5 


(60) 


is  the  mass  flow  rate  of  the  iC^  lubricant 

Cp^  is  the  specific  heat  of  the  i1"*1  lubricant  at  Temperature  T^ 

For  the  outer  chamber  surface,  the  appropriate  fluid  properties  are  those 
of  air,  and  the  temperature  difference  is  taken  between  the  chamber  wall 
and  the  external  ambient.  For  external  circulation  of  air  over  the 
chamber,  forced  convection  relations  must  be  used,  for  which  standard 
formulations  can  be  obtained  from  Kays^^. 

The  shaft  surface  coefficients  can  be  obtained  from  equations 
(52)  through  (58),  using  equation  (58)  rather  than  equation  (56),  if 
care  is  taken  to  use  mass  averaged  fluid  properties.  The  resulting 
prediction  will  be  higher  than  the  actual  cosfficient  value.  The  clear- 
ance  required  here  is  the  difference  between  shaft  and  chamber  radii. 

This  result  should  be  checked  against  the  prediction  of  equation (59) , 
which  provides  a  lower  bound  on  the  coefficient  value,  and  engineering 
judgement  used  to  choose  the  proper  valve  of  the  coefficient. 


The  heat  transfer  coefficient  for  the  external  surface  of  the 
loading  yoke  depends  on  the  geometry  of  the  member.  Generally,  this  can 
be  treated  as  a  fin  heated  at  the  root,  and  standard  correlations  applied, 
such  as  those  of  reference  [20]. 


Sample  Thermal  Analysis  Results 

The  results  of  the  thermal  analysis  of  a  conventional  inner 

g 

land  riding  roller  bearing  operating  at  2  x  10  DN  are  summarized  in  Table  3. 
The  guiding  shoulders  for  the  rollers  were  on  the  outer  race.  Two  cases 
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are  compared,  for  conditions  with  and  without  under  race  cooling. 

This  analysis  was  accomplished  using  specified  operating  values  of  IDC 
and  land  clearances,  and  hence  the  results  are  quite  comparable  for  cage 
slips  and  power  dissipations.  The  original  conditions  necessary  to  achieve 
the  specified  operating  conditions  were  computed  from  the  results  of  the 
thermal  analysis  and  considering  dynamic  influences,  and  are  summarized 
in  Table  3  along  with  results  pertinent  to  the  thermal  analysis. 

The  results  shown  in  Table  3  are  consistent,  showing  generally 
lower  temperatures  for  the  case  with  additional  cooling,  and  a  higher 
oil  temperature  rise  for  the  reduced  cooling  case.  Also,  with  under  race 
cooling  the  temperature  difference  from  the  inner  to  the  outer  race  is 
13.7°F,  as  compared  to  a  temperature  difference  of  34.7°F  for  the  case 
of  no  under  race  cooling.  These  results  are  in  general  agreement  with 
measured  values  for  comparable  bearing  operating  conditions. 

DISCUSSION  OF  APPLICATION  OF  EXPERIMENTAL  TRACTION  DATA  AND  ROLLER 
KINEMATICS 

Typical  traction  values  for  several  roller  bearings  from  65mm 
bore  to  140mm  bore  and  for  shaft  speeds  from  1.0  X  10^  DN  through  4.0  X  10^ 
DN  with  MIL-L-7808  oil  are  presented  in  Table  4,  for  unloaded  and  the 
maximum  loaded  roller  contacts.  All  tabulated  values  were  computed  with 
the  full  EHD,  fluid  mechanics,  dimensional,  and  thermal  models  described 
in  this  paper. 

A  roller  always  rotates  near  a  pure  rolling  condition  (low 
sliding  velocity)  at  the  outer  race  contact  of  a  roller  bearing  with  inner 
race  or  shaft  rotation.  Sliding  velocities  from  near  zero  to  3  in/sec  are 
typical  as  shown  in  Tables  4a  and  4b  for  outer-race  contact. 

An  unloaded  roller.  Table  4a,  rotates  about  its  own  axis  at  an 
angular  speed  slightly  less  than  ths  angular  speed  required  for  pure 
rolling  at  the  outer  race  contact.  Positive  sliding  (outer  race  relative 
surface  velocity  is  greater  than  the  roller  surface  velocity)  always 
occurs  and  is  in  the  range  of  near  zero  to  3  in/sec.  The  rolling  portion 
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of.  the  total  traction  value  is  significant  and  cannot  be  Ignored.  Indeed 
the  rolling  portion  can  even  predominate  as  shown  by  the  ratio  of  rolling 
traction  to  total  traction  in  Table  4a. 

A  loaded  roller  which  produces  a  cage  driving  force  rotates 
about  its  own  axis  at  an  angular  speed  slightly  greater  than  the  angular 
speed  required  for  pure  rolling  at  the  outer  race  contact.  Negative 
sliding  occurs  and  is  in  the  range  of  near  zero  to  -3  in/sec  as  shown  in 
Table  4b.  The  maximum  loaded  roller  and  most  of  the  loaded  rollers  in  the 
bearing  have  negative  sliding  at  the  outer  race.  Occasionally  a  lightly 
loaded  roller  or  rollers  near  the  entry  and  exit  from  the  radial  load  zone 
will  have  positive  sliding  at  the  outer  race  contact.  Such  rollers  result 
in  a  net  drag  on  the  cage  similar  to  completely  unloaded  rollers.  In  all 
cases  the  rolling  portion  of  the  total  traction  is  significant  and  cannot 
be  ignored. 

A  loaded  roller  always  has  positive  sliding  at  the  inner  race 
contact.  Typical  values  are  shown  in  Table  4c.  A  bearing  operating  at 
essentially  epicyclic  cage  speed  (no  cage  slip)  has  very  low  values  of 
sliding  velocity  (<  3  in/sec)  at  the  inner  race  contact  and  the  rolling 
portion  of  the  total  traction  is  significant.  A  bearing  operating  with 
measureable  cage  slip  can  develop  high  values  of  positive  sliding  velocity 
at  the  inner  race  contact  as  shown  in  Table  4c.  In  these  cases  the  sliding 
portion  of  the  traction  reaches  a  ceiling  value  and  again  the  rolling 
portion  of  the  total  traction  value  is  significant  and  cannon  be  ignored. 

A  high  speed  roller  bearing  with  elastohydrodynamic  oil  films 
cannot  operate  under  a  pure  rolling  condition.  An  examination  of  the 
traction  formulation  of  Eq.  (3)  shows  that  for  zero  sliding  vehocities  the 
traction  at  both  inner  and  outer  races  would  be  negative  and  no  equilibrium 
solution  can  occur  as  shown  in  Fig.  2.  Some  cage  slip  (however  small) 
must  occur. 
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RECOMMENDATIONS 


1,  Experimental  data  defining  the  rolling  portion  of  the  total 
traction  force  in  an  EHD  contact  is  needed.  The  accuracy  and  usefulness 
of  a  complete  systems  analysis  of  a  high  speed  roller  bearing  is  very 
dependent  upon  good  rolling  friction  data. 

Archfrd  and  Crook ^  have  stated  that  for  loads  which  develop 
essentially  parallel  surfaces  in  the  Hertzian  contact  and  at  speeds  of 
sliding  above  a  few  centimeters  per  second,  the  sliding  component  of 
traction  predominates  and  the  rolling  component  can  be  ignored.  This 
assumption  has  been  generally  accepted  and  has  led  to  an  emphasis  upon 
obtaining  twin  disk  traction  data  only  for  sliding  conditions  by  current 
researchers. 

Typical  twin  disk  machines  (with  equal  size  disks)  measure  only 
the  sliding  contribution  to  the  total  traction  in  an  EHD  contact.  Be¬ 
cause  of  disk  symmetry  the  data  is  also  symmetric  about  zero  sliding 
speeds,  i.e.,  the  traction  value  vanishes  at  zero  sliding  and  reversal 
of  tbv.  relative,  sliding  direction  reverses  the  direction  of  the  traction. 
The  rolling  portion  of  tue  total  traction  is  very  small  compared  to  max¬ 
imum  traction  values.  The  instrumentation  of  existing  traction  machines 
'is  not  designed  to  measure  small  magnitudes  of  rolll.  sig  traction,  indeed 
the  contribution  of  rolling  traction  is  calibrated-out  along  with  disk 
support  bearing  friction.  Symmetry  of  test  machine  geometry  leads  to 
symmetry  of  results.  Special  attention  to  non-symmetrical  bench  rig 
design  and  instrumentation  is  required. 

2.  Cage  dynamics  need  to  be  explored  and  better  defined  both 
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experimentally  and  analytically.  Bi-stable  operation  is  shown 
to  be  possible  by  the  present  analysis. 
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Figure  4.  Solution  Techniques  Flow  Chart 
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Table  2 

DESCRIPTION  OF  THERMAL  ANALYSIS  BY  NODES 
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TYPICAL  THERMAL  ANSLYSIS 


Bearing  t  Lubricant  Data 

Fitch  diameter  -  6.312  In. 

Roller  diameter  -  0.3  In. 

Operating  IOC  -  .002  In. 

Operating  radial  land  clearance  -  0.015  in. 

Lubricant  -  HIL-L-7808D 

Lubricant  flou  rate*  -  6  lb/sin  through  hearing  (both  runs) 
6  lb/ain  unc>r  race  (2nd  run  only) 
Lubricant  inlet  temperature  -  250*F 
Radial  load  -  333  lb 
Shaft  speed  -  14,2ft)  rpa 
011/air  ratio  In  bearing  -  152 

Temperature  Distribution  (Ref.  Pig.  9) 

Temperature,  *F 


No  Under-Race 

Under  Race 

Nodfc 

Location 

Coollns 

Coollna 

i 

Inner  shaft  coolant 

— 

— 

2 

Inner  shaft  surface 

361.5 

295.8 

3 

Outer  shaft  surface  at  bearing 

y  \  7 

293.6 

4 

Inner  race  at  shaft 

8 

310.7 

5 

Under-cut  surface  o»  inner  race 

362.5 

314.6 

6 

Under  race  coolant 

— 

265.9 

7 

Inner  race  beyond  roller  contact 

361.7 

317.1 

8 

Cage  at  the  land 

380.8 

370.3 

9 

C ««e  Inner  surface 

393.1 

383.1 

10 

Cage  outer  surface 

392.9 

383.0 

11 

InteriAl  lubricant  to  bearing 

330.7 

318.3 

12 

Rollers 

351.3 

339.2 

13 

Inner  race  at  rollers 

367.5 

324.8 

14 

Inner  surface  of  outer  race 

329.9 

30.-,  2 

15 

Inner  cooled  surface  cuter  race 

329.9 

307,: 

16 

Outer  race  lubricant 

— 

)  7 

Outer  cooled  surface  outer  race 

329.8 

306.2 

18 

Bearin*  cuter  diameter 

329.4 

TC2.6 

19 

Chamber  well  (inner  olde) 

309.0 

273,7 

20 

Loading  bolt  (strut)  root 

281.8 

.-60.1 

21 

Shaft  (external  to  bearing) 

335.3 

293.1 

No  Under-Race 

With  Under-Race 

Coollns 

Cooling 

Percent  cage  slip 

32.5 

33.7 

Pov-  r 

dissipation,  H 

P 

7.09 

7.15 

Initial  ehrinV  fit,  in. 

0.0082 

0  0088 

Machined  IOC,  m. 

0.0097 

0.0089 

Mounted  IOC,  in. 

0.0024 

0.0011 

Mithlncd  diametral  lend  clearance,  in. 

0.0037 

C.0035 

Mounted  diameter  lend  cleat  -;e,  in. 

0.0029 

0.0027 

Mired  lubricant  temperature  tt  steep,  *p 

330.7 
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Table  4 

TYPICAL  EHD  FILM  THICKNESS  AND  TRACTION  VALUES 


4a 

OUTER  RACE  CONTACT 
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•.Ml 
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i.MI 
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M 

14.M 
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•JOS 

Sill 

S.9M 

SM.fl 
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*  i.iu 

•  ? 

H 
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a.osi 

•.Ml 

-t.m 

4.1 

4.M4 

aus 

t.ili 

SM.il 

i.MM 

i.M* 

•  4.SS1 

OF  UNLOADED  ROLLERS 


t.Mt  Il.t  IM  IH-IM  tl  '  «.UM  tan  HI  ).• 

»•>»  #•'  IM  MtnlnMM  U1  t.Mt  l— •  »»  1.11 

IM  Ml-L-KM  m  I1M  tan  ■  1U  ?.» 

>■  «a-in«M  ut  t.ur  tan  im  t.n 

1.1M  (  ■  IX  IM-IM  IM  t.l.ta  itall  Ml  (.11 

t.m  t  <■  IM  HM-HH  IM  O.IM  tau  IM  «. n 

«.«t  m  MU  IM  t.lM  tan  Ml  e.ll 


4b 

OUTER  RACE  CONTACT  OF  MAXIHUM  LOADED  ROLLER 


f 

• 

f 

sum** 

Nu» 

^'(T)  » 

W.Um 

«ftu» 

1*4*1 

TMtUl 

its 

ru* 

1*4* 

Mi 

UUt 

•all 

MUS 

Imu 

Sd 

C4M 

|2hi.  telm. 

MM*. 

'.1  *  .1  *’ 

?4f 

gin 

ill  Im 

Mm 

Ms. 

««- 

•444*4 

im 

UM 

ImiumI 

s.i 

si.i 

SMI 

-1.14 

lJi.e 

4.M4S 

-  4.414 

-  i.MI 

*  s.m 

M.l 

ISO 

■U«Mt 

4» 

1.044 

MMr 

MLk— . 

MS 

.11  — 

t.ss 

S.M 

t.tsi 

MSI 

-4.M 

ssi.i 

OUf4 

•  .M4 

-  4.144 

•  4-JM 

St.t 

ISi 

WUL-S4M 

11! 

i.SM 

Iwi 

MIS 

i.M 

1.4 

•■•44 

ISM 

-i.44 

141.4 

4.414 

-  1.SC 

•  8»1M 

-  I.4IS 

4.1 

IS* 

IW-M 

14t 

i.SS4 

iMIt 

SSI 

i.iS 

>.c 

Si.M 

SPSS 

-s.m 

so.s 

4.4*4* 

•  i.lts 

-•.Ml 

•  S.MI 

ll.i 

ISi 

MS-4- TIM 

IM 

i.SM 

Umr 

III 

i.Si 

4.4 

•.•Mi 

Sill 

’•■•MS 

lll.S 

•■MM 

•  j.m 

•  i.ssi 

•  1. 444 

M 

JS4 

•U-WOt 

S«t 

f.M 

i4M> 

US 

i.ll 

4* 

M 

IMS 

^.4*i4'* 

SM.4 

W4 

-  .Ml 

-i.sn 

*•.»» 

i.i 

IS* 
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SM 

i.M 

Umi 

so 

i.M 

4,t 

4.M4 
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-4.4«1*** 

i.MM 

-  S.MI 

-i.tss 

-  I.U5 

s.i 

SM 

MU 

1*4 

f.UJ 

Hmt 

US 

e.u 

4c 

INNER  RACE  CONTACT  OF  MAXIMUM  LOADED  ROLLER 


Nnta  t*ntm 


ti 

SM 

SMI 

444.4 

MM 

i  ISM 

I.MI 

•  i.MI 

I.lM 

w.O 

ssf 

MUir**Mi 

ts 

4.0*4 

M4U 

OS 
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•  Ml 

SMS 

VIM 
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MM 

-•.MS 
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1S.1 
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MM 

MU 

•I 

IIS.I 
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14 
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S 
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IM 

I.M 
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*,» 

4  1 
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4SS4 
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SS4.i 

S.MI* 

MM 
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s.m 

4.S 

IM 

VA 

I.M 

taut 

*SS 

ta  V4 
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M 

MM 

:s» 

S*»4 

•.»** 

I.44S 

-  ».«’» 

*.}«£ 

•«.4 

SM 

atU'MI 

SM 

i.M 

iMt 

S*. 

6  S! 

•4 

•Ml 

ISO 

•  M 

tM.f 

4.414 

I.IU 

'  •  »*• 

I.SM 

*  * 

MS' 

***•*■»/ 

tilt 

S*4i 

IM 

i.M 

S*w» 

IS) 

4  S» 
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nomenclature 


Load  &  Speed  Equilibrium 


b 

Cr 

CF* 

d 

E 

"ym 

E* 

f(Te) 

Fddd 

F 


G 

G1,G2,G3 


,  m\l 

Phzh 


B1  "e v 
8K_ 


a  P, 


hz 


semi-width  of  Hertzian  contact  -  in. 

hollow  roller  ring  deflection  constant  -  in/lb. 

centrifugal  force  acting  on  roller  -  lb. 

roller  diameter  -  in. 

pitch  diameter  of  bearing  -in. 

2 

Modulus  of  elasticity  -  lb. /in 

2 

reduced  Modulus  *  E  /(l-o  ) 
ym 

coefficient  of  traction  at  temperature  Te 

contact  force  on  roller  from  guiding  race  shoulder  (lb.) 

cage  to  roller  force  -  lb. 

radial  external  load  -  lb. 

material  parameter  ■  aE’ 

parameters  to  determine  f(Te) 

(Shear  Rate  Effects) 


(Thermal  heating  effects) 
(Pressure  viscosity  effect) 
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lubricant  average  film  thickness  -  in. 

roller  polar  moment  of  inertia  lb-in-sec2 

cage  polar  moment  of  inertia  lb-in-sec2 

bending  moment  of  inertia  for  hollow  roller  -  in4 

roller  to  race  deflection  constant 

thermal  conductivity  of  lubricant,  BTU/o_-hr.-ft. 

F 

effective  length  of  roller  contact  -  in. 

roller  mass  -  lb-sec2/in. 

cage  fluid  drag  moment  -  lb-in. 

roller  fluid  drag  moment  -  lb-in. 

roller  load  -  lb 

order  of  iteration 

maximum  Hertz  contact  stress  -  lb. /in2 

diametral  clearancn  -  in. 

hydrodynamic  pressure  force  -  lb. 

dimensionless  form  of  Q  s  Q/(£E*R) 

dimensionless  heating  parameter ■  2n  U2/(K  T  ) 

e  f  e 

mean  radius  of  hollow  roller  -  in. 

radius  of  an  equivalent  cylinder  and  plane  -  in. 

thickness  of  hollow  roller  -  in. 

traction  in  an  EHD  contact  -  lb. 

temperature  of  lubricant  at  entrance  to  EHD  contact  -  °F 
fluid  entrainment  velocity  -  in. /sec. 
dimensionless  form  of  U  «  n  V/(E’R) 


I 

I 

I 

f 

% 

I 

I 

1 

J 

■f* 

I 
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TOTAL 


■  sliding  velocity  -  in. /sec. 

*  dimensionless  form  of  V  *  n  V/(E'R) 

■  dimensionless  load  ■  P/(iE'R) 

■  radial  displacement  of  inner  race  -  in. 

*  deviation  or  races  from  true  circular  form  -  in. 

*  total  number  of  rollers 


2 

■  pressure  viscosity  exponent  -  in  /lb. 

*  temperature  viscosity  coefficient  -  °F 

“  ratio  of  roller  diameter  to  bearing  pitch  diameter 

■  cogular  location  of  roller,  deg. 

*  thermal  reduction  factor  for  EHD  film  thickness 

*  angular  velocity  of  inner  race,  rad/sec. 

-  angular  velocity  of  cage,  rad/sec 

«  angular  velocity  of  roller,  rad/sec. 

2 

-  lubricant  absolute  viscosity  -  lb. -sec/in. 

■  coulomb  coeff  of  sliding  friction 

»  roller  -  outer  race  deflection  -  in. 

-  roller  -  inner  race  deflection  ~  in. 

,  * 

*  hollow  roller  deflection  from  CF  -  in. 

*  roller  diametral  expansion  from  -  in. 

*  total  deflection  of  an  unloaded  roller  -  in. 


relative  movement  of  races  at 


poison’s  ratio 


-  in. 
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t_  w 


u 

e 


SUBSCRIPTS 

*  inner  race  contact 

■  outer  race  contact 

*  roller  number  or  location 

■  an  unloaded  roller 

*  entrance  region  of  EHD  contact 


Fluid  Mechanics 


G 

Vji 


Vjk 

Ax 

C 

c’ 

CN 

FjG 

f 

fL 

L 


Surface  Area  in. 

2 

Total  contact  area  between  roller  and  guiding  shoulders  in. 
Area  of  icn  horizontal  lamina  in. 

Roller  and  cage  surface  area  within  typical  crossflow  passage 
exposed  to  crossflow  in. 

th  2 

Area  of  k  vertical  lamina  in. 

2 

Cross  sectional  area  of  a  typical  crossflow  passage  in. 

Radial  clearance,  in. 

Clearance  of  cage  pocket  and  roller,  in. 

Laminar  or  turbulent  correlation  factor  for  disk  rotating  in 
a  viscous  fluid 

Contact  force  between  roller  and  race  guiding  shoulder  lb 
Friction  factor 
Laminar  friction  factor 
2°aring  length 

Peripheral  length  of  a  typical  crossflow  passage  in. 

Moment  about  center  line  of  a  disk  rotating  in  a  viscous 
fluid  in-lb. 


-56- 

170 


1 

t 


N 


N, 


N, 


Re 

Ta 

r 


rHji 

ri 

r 

o 

rVjk 

S 

6 

P 

T 

Tadd 


Number  of  horizontal  lanina 
Reynold's  number  (Defined  in  Text) 

Taylor  number  (Defi 'ed  in  Text) 

Reference  radius  in. 

Distance  from  roller  center  to  horizontal  lamina  in. 

Inner  radius 
Outer  radius 

tH 

Distance  from  radial  ray  through  roller  center  to  k  lamina  in. 
Number  of  vertical  lamina 

2 

Crossflow  differential  pressure  lb/in. 

Drag  torque  in. -lb 

Moment  contribution  of  contact  force  to  roller  drag  moment 
in. -lb 


T 

eland 

T 

cside 

T 

rend 

T  , 
rcyl 


Fluid  drag  moment  on  cage  lands  in. -lb 

Fluid  drag  moment  on  cage  sides  in. -lb 

Moment  contribution  of  roller  ends  to  roller  drag  moment 
in. -lb 

Moment  contribution  of  roller  surface  to  roller  drag  moment 
in. -lb 


U  Mass  average  velocity  of  fluid  in. /sec 

V  Peripheral  crossflow  surface  average  velocity  in. /sec 

Velocity  of  roller  at  itl1  horizontal  lamina  in. /sec 

VRi  Velocity  of  race  at  ic^  horizontal  lamina  in. /‘sec 

^Vjk  Velocity  of  roller  at  k"*1  vertical  lamina  in. /sec 

uc  Coulomb  coefficient  of  friction  between  roller  and  race  guid¬ 
ing  shoulder 

2 

v  Fluid  Kinemetic  visosity  in.  /sec 

2  4 

p  Fluid  mass  density  lb  sec  /in. 

n 

Wall  shear  stress,  lb /in.* 
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Thermal  Analysis  and  Dimensional  Model 


Bk 

b 

C 


'Pi 


m,n 


ij 


IDC 

k 

kf 

9. 

s 

n. 


Ta 


corn! 


'Ta 


cr 


ava 


roller  inner  radius,  in. 

energy  rate  constants  in  formulation  of  thermal  problem, 
ref.  equation  (48),  Btu/sec. 

roller  outer  radius,  in. 

radial  clearance  between  cylindrical  surface  -  in. 

specific  heat  of  lubricant  or  coolent  associated  with  i 
thermal  node  -  B?.u/lb°F. 

Radiation  view  factor  between  m1"11  and  nth  surface 

gravitational  acceleration,  386.4  in. /sac2 

energy  transfer  between  nodes  i  and  j  Btu/sec 

convective  conductance  associated  with  ktk  surface  Btu/sec  - 
in.  2  °F 

bearing  internal  diametral  clearance  -  in. 

thermal  conductivity  -  Btu/sec  in.  °F 

thermal  conductivity  of  fluid  -  Btu/sec  in.  °F 

generic  length  of  cylindrical  surface  -  in. 

lubricant  or  coolent  flow  rate  of  i^k  thermal  node  —  lb/sec. 

total  number  of  component  dissipations,  external  energy  inputs, 
energy  transfer  to  coolent  or  lubricant  flows,  and  losses  at 
ith  thermal  node. 

total  subset  of  not  expiicity  containing  temperature. 

Taylor  Number  for  pure  conduction  through  a  fluid. 

critical  Taylor  number  at  speed  ratio  n. 

radius  of  inner  race  at  land  for  an  inner  land  riding  cage  or 
radius  of  outer  race  at  land  for  an  outer  land  riding  cage  - 
in. 

mixed  lubricant  temperature  -  °F 
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ref 


3 

Y 

AT 

At 

6 


ij 


T  . 
1 


-i 

Sb 


P 

n 

v 

c- 


temperature  of  i1"*1  thermal  node  -  °F. 

initial  uniform  reference  temperature  of  all  bearing  elements 

coefficient  of  thermal  expansion  -  in. /in.  °F. 

o 

thermal  diffusivity  -  in.  /sec. 

2 

kinematic  viscosity  of  fluid  -  in.  /sec. 

coefficient  of  volumetric  expansion  -  l/^F. 

Temperature  difference  between  i  and  j  thermal  nodes, 
W  -  °F. 

absolute  temperature  difference  between  surface  and  fluid. 

thermal  deformation  of  the  iC^  surface  of  a  bearing  member, 
(see  Figure  1)  -  in. 

til 

energy  error  function  at  ic  temperature  node  for  Newton- 
Raphson  solution  of  temperature  field  in  bearing  -  Btu/sec. 

til 

eraissivity  of  a  surface. 

coulomb  coefficient  of  friction  between  rollers  and  cage, 
ratio  of  cage  to  shaft  angular  relocities. 

Poisson's  ratio 

Boltzmann  radiation  constant  -  Btu/se.'  in.  “R" 
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